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DIMERS ON RAIL YARD GRAPHS 


CEDRIC BOUTILLIER, JEREMIE BOUTTIER, GUILLAUME CHAPUY, 
SYLVIE CORTEEL, AND SANJAY RAMASSAMY 


Abstract. We introduce a general model of dimer coverings of certain plane 
bipartite graphs, which we call rail yard graphs (RYG). The transfer matrices 
used to compute the partition function are shown to be isomorphic to certain 
operators arising in the so-called boson-fermion correspondence. This allows 
to reformulate the RYG dimer model as a Schur process, i.e. as a random 
sequence of integer partitions subject to some interlacing conditions. 

Beyond the computation of the partition function, we provide an explicit 
expression for all correlation functions or, equivalently, for the inverse Kaste- 
leyn matrix of the RYG dimer model. This expression, which is amenable to 
asymptotic analysis, follows from an exact combinatorial description of the op¬ 
erators localizing dimers in the transfer-matrix formalism, and then a suitable 
application of Wick’s theorem. 

Plane partitions, domino tilings of the Aztec diamond, pyramid partitions, 
and steep tilings arise as particular cases of the RYG dimer model. For the 
Aztec diamond, we provide new derivations of the edge-probability generating 
function, of the biased creation rate, of the inverse Kasteleyn matrix and of 
the arctic circle theorem. 


1. Introduction 

The two-dimensional dimer model is arguably the most studied exactly solvable 
model in statistical mechanics (note that it encompasses, in a sense, the equally 
well-known two-dimensional Ising model), see for instance |Kas67l Chapter 5] for a 
review of the seminal works of Kasteleyn, Temperley and Fisher, and the introduc¬ 
tion of [A1114j for a nice survey of the more recent literature. Dimer configurations 
are also known as perfect matchings in combinatorics and theoretical computer sci¬ 
ence. Actually, perhaps the oldest exact solution of a 2D dimer model is MacMa- 
hon’s enumeration of plane partitions [Mac04j . as these were later identified with 
lozenge tilings, or alternatively dimer configurations on the hexagonal lattice. 

Kasteleyn’s method allows to reduce the problem of computing the partition 
function (and the correlation functions) of the dimer model on any finite weighted 
planar graph (assuming that the dimers interact only through their hard-core re¬ 
pulsion) to the evaluation of a determinant (or Pfaffian) whose size is linear in the 
number of vertices of the graph. Under the usual assumption that the graph is 
periodic in two directions, one can then evaluate this determinant and take the 
thermodynamic limit to obtain the free energy, study phase transitions, etc. In 
this paper, we consider a dimer model on a new family of graphs, called rail yard 
graphs, which are periodic in one direction but not in the other. 
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One of our motivations is that the rail yard graph dimer model encompasses both 
the plane partitions mentioned above and another celebrated model, namely domino 
tilings of the Aztec diamond |EKLP92al IEKLP92b] (corresponding to, roughly 
speaking, dimer conhgurations on the portion of the square lattice htting into a 
large square tilted by 45°). What relates these two models is that they can be seen 
as Schur processes [OR 03] . that is to say random sequences of integer partitions 
whose transition probabilities are given by Schur functions. If the relation between 
plane partitions and Schur processes was explicited by Okounkov and Reshetikhin, 
the case of the Aztec diamond appears implicitly in |Joh02| and has, to the best 
of our knowledge, remained in such implicit form until |BCC14| , of which this pa¬ 
per is a continuation (see below). The interest of making the connection between 
dimer models and Schur processes explicit is that it allows to use an operator for¬ 
malism coming from the boson-fermion correspondence (see the references given at 
the beginning of Section]^ which is both powerful and intuitive, as the operators 
are nothing but transfer matrices or observables satisfying some particularly simple 
commutation relations. Furthermore it allows us to say that the RYG dimer model 
forms another situation, besides the 2D Ising model |Duh11| . where “bosonization” 
works at an exact discrete level. The rail yard graph dimer model corresponds 
essentially to the most general Schur process with nonnegative transition probabil¬ 
ities. 

Before describing our work in more detail, let us further discuss some history and 
background behind it. Bender and Knuth |BK72| made the link between plane par¬ 
titions and the Robinson-Schensted-Knuth correspondence, see also |Sta99l Chap¬ 
ter 7]. Okounkov jOkoOll rOko02j used the boson-fermion correspondence to de¬ 
fine the so-called Schur measure over integer partitions, and study its correlation 
functions. The Schur process |OR03l IOR07| is a time-dependent version of this 
measure, that can also be viewed as a system of particles with certain dynam¬ 
ics. It contains as a special case a generalization of plane partitions, namely plane 
partitions with an evolving “back wall”. Numerous papers followed on this sub¬ 
ject [Bor071IBF141 IBorlll IBMRT121 lBBB~*~14l IBFlSj and on its extension to the 
Hall-Littlewood and Macdonald cases [Vul07lIFW091 IVul09llOkalOlICSVlllIBC14| . 

In [BCCI4| . three authors of the present paper introduced a general class of 
domino tilings called steep tilings, encompassing both tilings of the Aztec diamond 
and the so-called pyramid partitions [Ken051 IYou09j . It was shown in [BCC14j that 
steep tilings also correspond to Schur processes and, using the vertex operator for¬ 
malism, their partition functions (of the “hook formula” type) were computed for a 
variety of boundary conditions. Since both (generalized) plane partitions and steep 
tilings are special instances of the Schur process, it is then natural to ask if there 
is a more general model of tilings or dimer coverings that would reformulate the 
Schur process in full generality, at least when the number of underlying parameters 
is finite. Such a model was sketched in [BCC141 Section 7], that can be viewed as 
a preliminary attempt at what we reach in the present paper. 

The rail yard graphs (RYG) that we introduce here are infinite bipartite plane 
graphs, obtained by the “concatenation” of column-shaped elementary graphs, and 
come with a family of admissible dimer coverings. The RYG dimer model is then 
a probability measure over such coverings. The elementary graphs can be of four 
types that correspond to the four possible types of “atomic” transitions in the 
Schur process. For the special families of RYG that correspond to the special 
families of Schur processes considered in [OR,031 IBCG14| , we recover generalized 
plane partitions and steep tilings, respectively. As we hope will be apparent in 
this paper, RYG provide a nice and natural formulation of the Schur process in 
terms of dimers, in a well-adapted system of coordinates, much simpler than the 



















































DIMERS ON RAIL YARD GRAPHS 


3 


one from |BCC141 Section 7]. Having shown the correspondence between rail yard 
graphs and Schur process, we can then apply the same classical tools as in |OB,03] to 
get explicitly the partition functions in a nice (hook-type) product form. Even more, 
we can interpret these partition functions in terms of a combinatorial parameter 
related to the flip operation on coverings. 

Beyond the partition function, we compute all the dimer correlation functions, 
which requires the introduction of suitably defined observables (or constrained 
transfer matrices) that enable us to localize, in the algebraic setting, a given set of 
dimers. To prevent any confusion, let us note that the particle correlations com¬ 
puted in [PRO3) for the Schur process, when translated in terms of RYG, give only 
a special case of this result. Indeed, as we will see, there are three kinds of dimers 
in a RYG, and particles correspond to one of the three kinds (so in our setting the 
correlations results of |OR03| only describe correlations between dimers of the first 
kind). Once the observables are constructed, we use classical fermionic tools such 
as Wick’s formula to evaluate the correlation functions in an explicit determinan- 
tal form. We also make the connection with the general Kasteleyn theory: it is 
a general fact that correlations between dimers on plane bipartite graphs have a 
determinantal form, underlaid by an inverse of the so-called Kasteleyn matrix of 
the model. For RYG, we show that the determinantal form we obtain by our ap¬ 
proach indeed gives an inverse Kasteleyn matrix, as was remarked in [OR07) for the 
case of skew plane partitions, see also |BF14l Section 5]. Our approach generalizes 
both this case and that of the Aztec diamond (for the so-called weighting), 
treated previously in |GY14| by a very tricky and somehow mysterious calculation. 
As further applications concerning the Aztec diamond, we rederive the so-called 
edge-probability generating function and biased creation rate, and the arctic circle 
theorem using the general saddle-point techniques of |OR03] . 

We now present the structure of the paper. Section is devoted to the basic def¬ 


initions (rail yard graphs in Subsection 2.1 their dimer coverings in Subsection 2.2 
flips in Subsection |2.3[) and to the statement of our main results, namely the ex¬ 


pression for the partition function (Subsection 2.41 and for the dimer correlation 
functions (Subsection |2.5| ). Section [^introduces bosonic operators (Subsection |3.1[ ) 
that act as transfer matrices in the RYG dimer model (Subsection |3.2[ ), allowing 
to compute efficiently the partition function (Subsection 3.3). Section [^ consid¬ 
ers fermionic operators (Subsection |4.1[ ) that play the role of observables in the 
RYG dimer model (Subsection |4.2[). Rewriting the correlation functions in the 
“Heisenberg picture” (Subsection |4.3| ), we derive their expression in the form of 
a determinant (Subsection |4.4[ ), before making the connection with Kasteleyn’s 
theory (Subsection 4.5). Section [^ discusses the previously known cases: plane 
partitions and lozenge tilings (Subsection 5.1) and steep domino tilings (Subsec¬ 
tion 5.2). In Sectionj^we address the specific case of the Aztec diamond, for which 
we provide new derivations of the edge-probability generating function and biased 
creation rate (Subsection |6.1[ ), of the inverse Kasteleyn matrix (Subsection |6.2[ ) and 
of the arctic circle theorem (Subsection |6.3[ ). Goncluding remarks are gathered in 
Section [^ Some auxiliary material is given in the appendix: a combinatorial proof 
of the bosonic-fermionic commutation relations (Appendix [^ and a rederivation 
of Wick’s formula (Appendix [^ . 

After the completion of this work, Alexei Borodin and Senya Shlosman informed 


us (by private communication) that graphs similar to those described in Section 2.1 


were suggested to them by Richard Kenyon, providing a dimer interpretation of the 
a,0-paths considered in [BSlOj . 
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2. Basic definitions and main results 

2.1. Rail yard graphs. We start by defining the underlying graph of our dimer 
model. We fix two integers £,r such that £ < r, and denote by [£..r] the set of 
integers between £ and r. We then consider two binary sequences indexed by the 
elements of \£..r\-. 

• the LR sequence a = (ag, a^+i,..., a^) S {L, 

• the sign sequence b = {bi, 6^+i ,... ,bj.) G {+, — 

The rail yard graph associated with the integers £ and r, the LR sequence a and the 
sign sequence &, and denoted by RYG(£, r, a, b), is the bipartite plane graph defined 
as follows. Its vertex set is [2£ — 1..2r + 1] x (Z + 1/2), and we say that a vertex 
is even (resp. odd) if its abscissa is an even (resp. odd) integer. Each even vertex 
{2m, y), m G [£--r], is then incident to three edges: two horizontal edges connecting 
it to the odd vertices {2m — l,y) and (2m +1, y), and one diagonal edge connecting 
it to 

• the odd vertex (2m — 1, y + 1) if am = L and bm = +, 

• the odd vertex (2m — l,y — 1) if o™ = L and bm = 

• the odd vertex (2m + 1, y + 1) if Om = R and bm = +, 

• the odd vertex (2m + 1, y — 1) if am = R and bm = —■ 

Hopefully, this explains our motivations for using the symbols L, R and +, —. Draw¬ 
ing the edges straight, the graph is indeed bipartite and plane by construction. For 
e an edge, we write e = (a, j5) to mean that a is the even endpoint of e, and /3 
its odd endpoint. For v a vertex, we will denote by u’' S Z its abscissa, and by 
S Z -b 1/2 its ordinate. 



'■ Lj^ • R- • R^ • L- ■ 


2f- 1 


2ril 


O odd vertex 


■ even vertex 


Figure 1. The rail yard graph associated with the LR sequence 

RLRRLR and the sign sequence -b H-1-(with .^ = 1, r = 6). 

It is infinite and periodic in the vertical direction, but finite in the 
horizontal direction. 


Figure displays the rail yard graph associated with the LR sequence RLRRLL 

and the sign sequence -b H-1-(with ^ = 1, r = 6). Observe that a rail 

yard graph is infinite and 1-periodic in the vertical direction. When £ = r, the LR 
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and sign sequences both consist of a single element, and the corresponding rail yard 
graph, which is said elementary, is of one of four possible types, see Figurej^ Given 
two rail yard graphs RYG(£, r, a, 6) and RYG(£', r', o', 6') such that £' = r + 1, we 
define their concatenation by taking the union of their vertex and edge sets. It is 
nothing but the rail yard graph RYGft?. r', aa'. bb'^ where aa! and W denote the 
concatenations of the LR and sign sequences. Glearly, a general rail yard graph is 
obtained by concatenating elementary ones. 



Figure 2. The four elementary rail yard graphs L+, L—, R+, 
and R— (up to horizontal translation) 


The left boundary (resp. right boundary) of a rail yard graph consists of all odd 
vertices with abscissa 2£ — 1 (resp. 2r + 1). Vertices which do not belong to the 
boundaries are said inner. When drawn in the plane, the graph delimits some faces, 
and the bounded ones are called inner faces. Note that inner faces may be incident 
to 4, 6 or 8 edges. Finally, observe that our definition works equally well if we take 
£ = —oo and/or r = +oo, thus considering infinite LR and sign sequences. In that 
case, the rail yard graph “fills” either the whole plane or a half-plane, boundaries 
being sent to infinity. 

2.2. Admissible and pure dimer coverings. We now turn to the characteriza¬ 
tion of the configurations of our dimer model. Given a rail yard graph RYG(£, r, a, b) 
with £, r finite, an admissible dimer covering is a partial matching of this graph such 
that: 

• each inner vertex is covered (i.e. matched), 

• there exists an integer N >0 such that: any left boundary vertex {2£—l, y) 
is covered for y > N and uncovered for y < —N, any right boundary vertex 
(2r + l,y) is covered for y < —N and uncovered for y > N, 

• only a finite number of diagonal edges are covered. 

A pure dimer covering is an admissible dimer covering for which the second property 
above holds for N = 0: in other words the uncovered vertices are precisely the left 
boundary vertices with negative ordinate and the right boundary vertices with 
positive ordinate (see Figure]^. The fundamental dimer covering is the pure dimer 
covering where no diagonal edge is covered (it is not difficult to check its existence 
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and uniqueness e.g. by induction on r — £). Observe that any admissible dimer 
covering coincides with the fundamental dimer covering outside a finite region. An 
elementary dimer covering is an admissible dimer covering of an elementary rail 
yard graph (see Figure]^. 




Figure 3. Some elementary dimer coverings. The underlying el¬ 
ementary rail yard graph has type R+ in the first two cases, and 
R—, L+, and L—, in the three others, from left to right. Outside 
of the displayed region, the configuration is identical towards the 
top (resp. bottom) on each horizontal level to what it is on the 
topmost (resp. bottommost) displayed level. 


Similarly to rail yard graphs, admissible dimer coverings behave nicely with 
respect to concatenation. More precisely, consider two rail yard graphs G = 
RYG(£, r, a, 6) and G' = RYG{£',r',a',b') which are concatenable (i.e. £' = r + 1) 
and let GG' be their concatenation. Let C and C' be admissible dimer coverings of 
respectively G and G'-. we say that G and G' are compatible if, for each y G Z-1-1/2, 
the vertex (2r + l,y) = (2£' — l,y) is covered in G if and only if it is not covered 
in C'. In that case, by taking the union of C and G', we obtain an admissible 
dimer covering of GG/ which we denote by GG'. Conversely, any admissible dimer 
covering can be decomposed as the concatenation of elementary dimer coverings 
which are sequentially compatible. 

It is also interesting to consider the limiting cases £ = —oo and/or r = -boo, 
which requires a slight adaptation of our definitions. An admissible (resp. a pure, 
resp. the fundamental) dimer covering is then a matching such that each inner 
vertex is covered, and such that there exists finite integers £',r' such that : 

• inside the strip [2£' — l,2r' -b 1] x K, we see an admissible (resp. a pure, 
resp. the fundamental) dimer covering in the previous sense, 

• outside this strip, all covered edges are horizontal. 

(Note that this definition works in all situations: it coincides with the previous one 
when £,r are both finite.) 

Our motivation for considering pure dimer coverings of rail yard graphs is that 
we recover several well-known dimer models as specializations. For instance, taking 
r — £ = 2n, a LR sequence of the form LRLRLR ■ ■ ■ and a sign sequence of the form 
H-1-1— • • •, the corresponding pure dimer configurations are in bijection with 
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, O '“O -“O '“O -“O , 

■ ■ L_ ■ i?_ ■ 

2e - 1 2r+l 


Figure 4. A pure dimer covering of the rail yard graph of Figure [l] 

domino tilings of the Aztec diamond of size n. We also recover plane partitions and 
so-called pyramid partitions, which requires taking i = —oo and r = -|-oo: plane 
partitions are obtained by taking a constant LR sequence and a sign sequence of 

the form ••• + + + -1- •••, while pyramid partitions are obtained by taking 

an alternating LR sequence (• • • LRLRLRLR • • •) and the same sign sequence. We 
will discuss these specializations in greater detail in Section 

2.3. Flips. We now define a local transformation on admissible coverings called 
the flip. Let G be a rail yard graph, C be an admissible covering of G, and let / 
be an inner face of G. If exactly half of the edges bordering / belong to G, then 
removing these edges from G and replacing them by the other edges bordering / 
gives another admissible covering C of G. The operation that replaces G by C is 
called the flip of the face /, see Figure 

Figure 5. Examples of flips. All these flips are positive when the 
transformation is made from left to right. 

We say that the flip of an inner face / is positive if after performing the flip, 
the edges of / that belong to the covering are oriented from odd to even vertices in 
counterclockwise direction around /. The flip is negative otherwise. For example, 
each flip displayed on Figure is positive when performed from left to right. It 
follows from |Pro021 Theorem 2] that the positive flip relation endows the set of all 
pure coverings of a given rail yard graph with a distributive lattice structure. In 
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particular, each rail yard graph has a unique minimal pure covering from which all 
other ones can be reached using positive flips only. The minimal covering is the 
only pure covering on which no negative flip is possible. Using this criterion one 
easily checks that the minimal covering coincides with the fundamental covering 
defined above. The flip distance between two coverings is the minimal number of 
flips needed to go from one to the other. When one of the two coverings is the 
fundamental one, the flip distance is realized by a sequence that uses positive flips 
only. 


2.4. Enumeration. Our main enumerative result is an expression for the partition 
function of the RYG dimer model, which we now define. Consider a rail yard graph 
G = RYG(£, r, a, b), and a sequence of formal variables x = {xi, x^+i,..., Xr), where 
possibly t = —oo or r = +oo. The weight of an admissible dimer covering C of G 
is then defined as 


( 1 ) wiC)=Ylxf^^^ 

i=£ 

where di{C) is the number of diagonal dimers in column i (i.e. the number of covered 
diagonal edges incident to an even vertex with abscissa 2i). This weight is well- 
defined since ^di{C) is finite by the definition of an admissible dimer covering. 
The partition function of the multivariate RYG dimer model, denoted Z{G;x), is 
then the sum of the weights of all pure dimer coverings of G. 


Theorem 1. The partition function of the multivariate RYG dimer model reads 


(2) Z{G-,3f)= n 

i<i<j<r 
bi—-\-,bj= — 


where 


(3) 


Zij 


1 XiXj if di ^ Clj j, 

(1 - ifai = aj. 


Remark 2. The partition function is always a well-defined power series in the xfs: 
indeed, all but finitely many factors contribute a factor 1 to the coefficient of a 
given monomial in (H). 

An interesting specialization is the q-RYG dimer model: given a formal variable 
q, we attach to each configuration a weight q'^ with d its flip distance to the funda¬ 
mental one. As explained in Section [3.3| below, this can be achieved by taking, for 
all i G [f..r], Xi = qRf hi = —, and Xi = Xjcf if = -I-, with q an indeterminate. 
A caveat is that, when £ or r is infinite, this specialization may be ill-defined since 
an infinite number of monomials in the xfs might specialize to the same monomial 
q‘^. A sufficient condition for the specialization to be well-defined is the following 
finiteness condition on the sign sequence: 

• if £ = — oo, then there exists £' finite such that hi = + for all i < £', 

• if r = -too, then there exists r' finite such that hi = — for all i > r'. 

(This condition is essentially necessary, because any initial run of — or final run 
of -I- in the sign sequence does not contribute to the partition function, and can 
be removed without loss of generality: any pure dimer covering coincides with the 
fundamental dimer covering in the corresponding regions.) 
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Theorem 3. Assuming that the finiteness condition holds, the partition function 
of the q-RYG dimer model is 

(4) Z{G-q)= n 

bi=+,bj=- 


(5) 


1 + q^ ^ ifoif^ Qj, 

ifai = aj. 


The proofs of Theorems and are given in Section ^ 


Remark 4. The product form Q is strongly reminiscent of a hook-length formula, 
upon interpreting the sign sequence as describing the shape of a (possibly infinite) 
Young diagram, see Figurej^ The finiteness condition ensures that there are finitely 
many “hooks” of a given length, and hence that Q is a well-defined formal power 
series in q. 



Figure 6. The interpretation of Theorem]^ as a hook-length for¬ 
mula. Displaying the sign sequence 6 as a lattice path whose hor¬ 
izontal (resp. vertical) steps correspond to -b (resp. to —), one 
obtains a Young diagram whose boxes are indexed by pairs (i,j) 
such that i < j, bi = +, bj = —. The quantity h = j — i is the 
“hook-length” of the box. In Q, each box gives rise to a multi¬ 
plicative factor Zij whose value l + q^ or (1 — q^)~^ is determined 
by comparing the two terms and Oj of the LR sequence a. 


2.5. Correlations. So far we have introduced the partition function of the RYG 
dimer model, which depends on a sequence of formal variables x in general and on 
a single variable q in the flip specialization. For a probabilistic or statistical physics 
interpretation, one shall rather consider the xfs or q as nonnegative real numbers 
such that the sum Z{G-,x) of the weights w{C) over all pure dimer coverings C is 
convergent. As apparent from Theorem this is the case if and only if 

(6) XiXj < 1 for all i < j such that Oi = aj, bi = + and bj = — 
and, when ^ or r is infinite, 

(7) E XiXj < oo. 

£<i<j<r 

bi=+.bj= — 

In the g-RYG dimer model, ^ is satisfied whenever <7 < 1, and Q amounts to the 
finiteness condition defined before. 

Assuming that the RYG dimer model is well-defined, that is to say Q and 0 
are satisfied, we may interpret w{C)/Z{G;x) as the probability of the pure dimer 
covering G. For a finite set E of edges, we denote by Pg-,x{E) the probability 
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that all the edges of E are covered by a dimer. Our main probabilistic result is 
an explicit determinantal expression for Pc-xiE), which requires to introduce some 
notations. For k, k' two integers, we set 


n {^+x^z) 


( 8 ) 


Fk{z) = 


2i<k 


2j>k 


I — Zl 
Z 


n (1 


n 


r.{ai,bi) = {L,+) 
2i<.k 


j:{aj,bj) = {R,-) 
2j>k 


1 + 

Z 


and 


(9) 


Gk,k'{z,w) 


Fk{z) ^/zw 
Fk' (w) z -w 


Note that all the products in ([^ are convergent by Q, hence Fk{z) is a meromor- 
phic function on the whole complex plane, whose all zeros and poles are on the real 
axis. 

For a, P two vertices of G such that is even and is odd, we set 


( 10 ) 


Ha,/? 


1 

(2z7r)^ 




Ga^^^^{z,w) 


dz dw 
z w 


where the contours must satisfy the following conditions: (i) Gz should encircle 
0 and all the negative poles of Fa^{z), but not the positive ones; (ii) G^ should 
encircle 0 and all the positive zeros of Fi 3 ^{w), but not the negative ones; (iii) Cz 
and Cuj should not intersect, and Gz should surround if and only if < P^. We 
shall check in Section 4.4 that the assumptions § and Q imply that such contours 
always exist but, at this stage, let us mention their intuitive interpretation: Ca,p 
is obtained by extracting the coefficient of z°‘^w~^^ in Go,^^py^{z,w), when we treat 
each factor (1 — Xiz)~^ as a power series in z, each factor (1 + Xjjz)~^ as a power 


series m z 


-1 


-1 


each factor [l+Xiw) ^ as a power series in w, each factor (1 — Xj/w) 
as a power series in and finally we expand Z^zwK^ — w) = 'ZkeN+i/ 2 ('^/ 
if < P^, or ^/zwl{z - w) ^ Efc 6 N+i/ 2 (w^/^)”'' if 

We are now ready to state our theorem. Recall that, for an edge e = {a, P), a 
and P are assumed to be respectively the even and the odd endpoint of e. 


Theorem 5 (Dimer correlations). Let E = {ei,..., e^} be a finite set of edges of 
RYG(.^, r, a, 6), with = {ai,Pi). Then, we have 


( 11 ) 


Pg-,x{E) = det (Co 


with H{E) the number of horizontal edges in E whose right endpoint is at an even 
abscissa, x— = xf’- ■ ■ ■ xf^ with Uk the number of diagonal edges in E in column k, 
and C defined as in ([I^ . 


The proof of Theoremj^is given in Section]^ where we also prove that the infinite 
matrix C, with rows indexed by even vertices a and columns by odd vertices P, is 
an inverse of the Kasteleyn matrix of the rail yard graph for a suitable Kasteleyn 
orientation (see Theorem 17). Applications will be discussed in Section]^ 


3. Bosonic operators 

The purpose of this section is to establish Theorems and This is done natu¬ 
rally by the transfer-matrix method, which here consists in decomposing the pure 
dimer coverings we want to enumerate into a sequence of compatible elementary 
dimer coverings. It turns out that the transfer matrices are isomorphic to certain 
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operators arising in the so-called boson-fermion correspondence. For more details 
on this latter subject, we refer the reader to one of the many references available 
in the mathematical physics literature, for instance |Kac901 Chapter 14], |MJD00j . 
[OkoOll Appendix A], |OR07| , |Tinll) and [AZ13| . 

We start by giving the necessary reminders in Section |3.H then make the con¬ 
nection with rail yard graphs in Section |3.2[ and finally complete the proofs of 
Theorems [T] and [3| in Section 13.31 


3.1. Reminders. An integer partition, or partition for short, is a nonincreasing 
sequence A = (Ai)i>i of integers which vanishes eventually. The size of a partition 
|A| =E*>i \i. We say that two partitions A and p are interlaced, and we write 
A ^ or p, ^ A, if we have 

(12) Ai > > A2 > /r2 > A3 > • • • 

In the well-known pictorial representation in terms of Young diagrams, this means 
that the skew shape X/p is a horizontal strip, see e.g. |Sta991 Chap. 7] for more 
precise definitions. To a partition A we may associate its conjugate X' , whose Young 
diagram is the image of that of A by a reflection along the main diagonal. In more 
explicit terms, we have A' = #{j > 0, Xj > i}. Note that the relation X >~ p 
amounts to 


(13) A' — /r' G {0,1} for all i >1. 

The bosonic Fock space, denoted B, is the infinite dimensional Hilbert space 
spanned by orthonormal basis vectors |A) where A runs over the set of integer 
partitions. Here we will use the bra-ket notation so that (A| denotes the dual 
basis vector. For x a formal or complex variable, we introduce the operators 
rL+(ai),ri_(a;),rij+(a;),rij_(a:) whose action on basis vectors reads 

Tl+{x)\X)= Y. x\^\-\>^\\p), Tr+{x)\X)= Y 

rr.-ix)\x)= Y r«-(*)|A)= E 

These operators are sometimes called (half-)vertex operators. Let us mention that, 
in the literature, rL±(x) is often denoted r±(a;), see e.g. |Oko01| . while ri^±(a;) is 
sometimes denoted Ffj. {x) [YoulOl IBCC14| . Observe that we have 

(15) ri+(x)|0) = r«+(x)|0) = 10), (0|ri_(a:) = (0|r^_(a;) = (0| 

where 0 denotes the empty partition. Note also that F^- (resp. F^-) is the dual 
of rL+ (resp. rfl_|_), and that r/j_|_ (resp. F^j-) is conjugated to ri_|_ (resp. F^-) 
via the involution lu oi B sending |A) to |A'). 


Remark 6. For 01,02 G {L,R}, the product rai-(a;i)ra 2 _|_(a; 2 ) is clearly well- 
defined, because its coefficient between two states (A| and \p) involves only a finite 
sum. The same is true for roj_|_(xi)r(j 2 _(x 2 ) when oi ^ 02 (observe that the “in¬ 
termediate” partitions cannot get too large). Infinite sums arise when considering 
raj_|_(a;i)ra 2 _(a; 2 ) with oi = 02 , but its coefficients are power series in xi and X 2 , 
which are convergent for |xia: 2 | < I as apparent from the following proposition. 


Proposition 7 (Commutation relations). For 01,02 G {L,R}, we have 

(1 - a;ia;2)“^ ra2_(o:2)Faj + (a;i) if ai = 02 , 


(16) 


rai+(a;i)ra2-(a;2) = 


( 1 - 


■ X1X2} 


ra2-(x2)rai + (a;i) f/01^02, 


while rai + (a:i) commutes with rtj 2 +(a; 2 ), and r£jj_(a:i) commutes with Ta 2 -{x 2 )- 
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Proof. See for instance [YoulOl Lemma 3.3] for an algebraic proof, and BBB+ldl 
Section 3] for a bijective proof of ([T^. The celebrated Bender-Knuth involution 
|BK721 pp. 46-47] implies that ra±(a;i) commutes with ra±(a; 2 ) for a = L or R. 
That r L±ixi) commutes with r/j±(a; 2 ) is also well-known, but for completeness let 
us here sketch a short proof: for two partitions A, /r, one sees easily that the two 
sets {v ■ ^ ^ V, v' ^ X'} and {v ^ v', v ^ X\ are nonempty if and only if 
A/^ is a skew shape containing no 2 x 2 square. In that case, both sets have the 
same cardinality 2*^, where C is the number of connected components of A//r, and 
one easily constructs a bijection between them proving the wanted commutation 
relation. Another byproduct of this bijection is that 


(17) 


rL±(x)ri^±(-a:) = ri^±(-a;)rL±(a;) = 1. 


□ 


Given two symbols • and o (called respectively black and white marbles), a Maya 
diagram |MJD00j is an element m of {*, 0 }^+^/^ such that m^, is eventually equal 
to • for k —>■ —00 and to o for /c —>■ -l-oo. It then not difficult to check that the 
quantity 

(18) c = #{fc > 0, mfe = •} - #{fc < 0, nifc = 0 } 

is a finite integer, and we call it the charge of m. Let A:i > fc 2 > • ■ ■ be the 
positions of • in m enumerated in decreasing order, and let Xi = ki — c + i — 1/2: 
it is easily seen that A is a partition and that the correspondence m i—> (A, c) is 
one-to-one, the pair (A, c) being called a charged partition. Observe that we may 
extend the involution uj to charged partitions by setting a;(A, c) = (A', —c), and this 
corresponds on Maya diagrams to performing a reflection across 0 and exchanging 
• and o: in other words, m is sent to m' such that {mfe,m'_ -J = {•> 0 } for all 
k G Z + 1/2. By a slight abuse, we still denote by 0 the Maya diagram of charge 0 
corresponding to the empty partition. 

The fermionic Fock space, denoted F, is the infinite dimensional Hilbert space 
spanned by orthonormal basis vectors jm) where m runs over the set of all Maya 
diagrams. For c G Z, let Fc C F denote the subspace spanned by Maya diagrams 
of charge c, so that F = (Bc&zFc. By the bijection between Maya diagrams and 
charged partitions, each Fc may be canonically identified with B. This defines 
the action of the bosonic operators Tand F ji± on F, leaving each subspace Fc 
invariant (by a slight abuse we keep the same notations for the operators acting on 
this larger space, and note that the commutations relations of Proposition[^remain 
valid). We now end this section devoted to reminders, leaving the discussion of 
fermionic operators to Section [4.1[ 


3.2. Interpretation as transfer matrices for RYGs. The purpose of this sec¬ 
tion is to explain how the operators F/F ii± may be identified with dimer transfer 
matrices. The key observation is that Maya diagrams describe the boundary states 
in our model. More precisely, let us consider an admissible dimer covering C of 
G = RYG(£, r, a, &). We define the left boundary state l(C') of C by setting, for all 
k gZ + 1/2, 


(19) 


KC)k 


o if (21! — 1, A:) is covered by a dimer, 
• otherwise. 


It is a Maya diagram by the definition of an admissible dimer covering. Similarly, 
the right boundary state r(C) of C is the Maya diagram defined by 


( 20 ) 


r{C)k 


• if (2i! -I- 1, A:) is covered by a dimer, 
o otherwise. 
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for all k G Z + 1/2. See Figure [^a). A pure dimer covering has both boundary 
states equal to 0. Note that, if G' is a rail yard graph which is concatenable after 
G, and if C' is an admissible dimer covering of G', then G and G' are compatible 
if and only if r(G) = l(G'). We may now state the main result of this section. 


Proposition 8 (Transfer matrix decomposition). Given a rail yard graph G = 
RYG(.^, r, a, 6) with i,r finite, and two Maya diagrams 1 and r, the sum of the 
weights Q of all admissible dimer coverings of G with left boundary state 1 and 
right boundary state r is given by 

(21) Z{G,\,v,x) = {\\Tatbt{xt)'^ai+ibt+A^t+i) ■ ■ ■ra,.h,.(a;^)|r). 

In particular, the partition function reads 


( 22 ) 


Z{G-,x) = {f!l\Ta,,bfixe)Tai^ibi+iixi+i) ■ ■ ■Ta,.b,.iXr)\9). 


Proof. Note that ( |22[ ) follows from 
considering pure dimer coverings. 


by taking 1 = r = 0, which amounts to 


We first verify (211 for I = r, i.e. when G is an elementary rail yard graph. Let 
us here treat the case = L, &^ = + (displayed third on Figurej^c)) and leave the 
other cases to the reader. Let Si < S 2 < • • ■ (resp. < ^2 < ■ ■ •) t>e the positions 
of o in 1 (resp. r) enumerated in increasing order. Then, we claim that both sides 


of ( |^ are equal to x 

(23) 

(24) 




if the two conditions 


Si — ti G {0,1} for alH > 1 


i>l 


hold, and that both sides vanish otherwise. Indeed, on the one hand, it is not 
difficult to check (see Figure]^ that there is at most one elementary dimer config¬ 
uration with prescribed boundary states 1 and r, and that there is exactly one such 
configuration (containing ^('5* ~ diagonal dimers) if and only if (23) and (24) 
hold. On the other hand, let A and p, be the integer partitions associated with the 
Maya diagrams 1 and r: the quantity (l|Fi+(a;^)|r) is equal to if the two 

conditions 


(25) 

(26) 


A ^ /r 

1 and r have the same charge c 


hold, and it vanishes otherwise. But we have A' = c -I- z — 1/2 — and = 
c -I- z — 1/2 — for all z > 1 hence, in view of (13), we find that the conditions 


(23)-(24) amount to (25)-(26l, and then that |^| — |A| = ~ ^i) wanted. 

It remains to verify (21) for £ < r, which may be easily done by induction: it 


suffices to observe that, if G' is the rail yard graph obtained by removing the last 
“strip” of G, then any admissible dimer covering G of G with boundary states 1, r 
is uniquely decomposed into a pair formed by an admissible dimer covering G' of 
G' with boundary states 1, m, for some Maya diagram m, and an elementary dimer 
covering E with boundary states m, r, such that w{G) = w{G')w{E). □ 


3.3. Proof of enumeration results and computation of the partition func¬ 
tion. We are now ready to prove Theorems and The first one is a direct 
consequence of the formalism developed above, whereas the second deserves an 
inspection of the different types of flips in rail yard graphs. 


Proof of Theorem\^ We simply have to evaluate the right-hand side of ( [2^ , which 
can be done as in in IOR071 Section 4.1], [YoulOl Section 4] or |BCC141 Section 
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Figure 7. (a) The rules ([T^-(20| that define the two Maya di¬ 
agrams l(C') and r(C'). (b) An admissible dimer covering of the 
elementary rail yard graph L+. White marbles are numbered as 
in the proof of Proposition (c) If the left boundary vertex at 
ordinate Si is covered by a horizontal (resp. diagonal) dimer, then 
the right boundary vertex at ordinate Si (resp. Si — 1) is necessarily 
uncovered. Therefore there is a white marble at this ordinate on 
the right boundary, and induction implies that it is the i-th white 
marble on the right boundary. This proves that Si — U G {0,1} 
and that ~ ^i) equal to the number of diagonal dimers. 

Moreover, since dimers incident to the left boundary can be recov¬ 
ered from the knowledge of the sequences (si) and (ti), it is clear 
that there is at most one dimer configuration with given boundary 
states 1 and r, and that there is one if and only if (23) and (24) 
hold, (d) Condition (231 can be interpreted by saying that white 
marbles are “jumping downwards” by 0 or 1. The cases of the 
elementary graphs L—, R+, and R— have similar interpretation, 
respectively with white marbles jumping upwards, black marbles 
jumping upwards, and black marbles jumping downwards, in all 
cases by 0 or 1. 


5.1]: first observe that, by (15), for any k,m, any (ci, C 2 ,... Cfe+m) G 
one has: 


fc+m 


(27) 


inrc.-(^.) n r.,+(z,)|0) = i, 




the Zi being formal variables. 

Now, by applying successively the commutation relations of Proposition one 
can transform (22) into a scalar product of this form, up to a multiplicative prefac¬ 
tor, by moving to the left all the operators Ta^bj (xj) such that bj = —. In order to 
do that, we have, for each £ < i < j < r such that — + and bj = —, to transform 
the product of operators ro^+(a:i)rf,^._ (xj) into the product Tb^-{xj)Ta- + {xi). For 
each such transformation, we obtain a multiplicative contribution given by (16), 
and the result follows. □ 
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Proof of Theorem^ As explained before the statement of Theorem we have to 
prove that, specializing for i G [Z..r] the weights Q to Xi = 9 ® if = — and 
Xi = l/g* if bi = + amounts to attaching to each conhguration a weight q‘^, where 
d is the flip distance to the fundamental configuration. 

First, this is true for the fundamental configuration that receives a weight 1 in 


both cases. Second, since by Propp’s theory (recalled in Section 2.3) each shortest 


path from the fundamental state to any configuration is realized using positive flips 
only, it is enough to check that, in this specialization, each positive flip increases 
the weight of a configuration by a factor of q. 




R+ 




Rr 


y 






"a . 

column i column i + 1 

(a) 


yy 7 

X TTo 

column i column i + 1 

ib) 


^ 7 
y 77 
77 7 
7 7X 

column i column i + 1 

ic) 


Figure 8 . (a) The 16 possible face types of rail yard graphs are 
obtained by matching one of the half-face types on the left with 
one on the right, (b-c) The dimer configuration around half-faces 
before (b) and after (c) a positive flip. 


Consider an inner face / in a rail yard graph. Then / is made by the union of 
two half-faces as shown on Figure [^a). Each of these two half-faces is incident to a 
diagonal edge, one in column i, and one in column f -|- 1, in the sense of Section|2.4 


for some i G [l..r — 1]. Then, a case inspection (see Figure |^b-c)) shows that the 
following is true: when performing a positive flip on /, the number of diagonal 
dimers on column i increases (resp. decreases) by 1 if 5^ = -|- (resp. bi = —), and 
the number of diagonal dimers on column * -|- 1 decreases (resp. increases) by 1 if 
bi+i = -I- (resp. = -). 

Therefore this flip multiplies the weight Q of the configuration by a factor of: 


(28) 


i / ® 


- 1-1 ’ 


where we identified -I-, — with -|-1, —1 respectively. But, in the specialization we are 

bi 


considering, we have xf = q for all j G [h.r], so (28) is equal to q and the proof 
is complete. □ 


4. Fermionic operators 
The purpose of this section is to establish TheoremWe start in Section 4.1 by 


recalling the definitions and basic properties of fermionic operators. In Section 4.2 


we show that these operators can be used to construct constrained transfer matri¬ 
ces, that is transfer matrices enumerating dimer configurations containing a given 
subset of edges. We rewrite the product of constrained transfer matrices in another 
convenient form in Section[4(^ and complete the proof of Theoremj^in Section 4.4 


Finally, we elucidate the connection with Kasteleyn’s theory in Section 4.5 
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4.1. Reminders. Recall that the fermionic Fock space T^ introduced at the end 
of Section |3.11 is the infinite dimensional Hilbert space spanned by orthonormal 
basis vectors |m) where m runs over the set of all Maya diagrams. For k G Z + 1/2, 
we define the fermionic operators ipk and tpl. (also called creation/annihilation 
operators) through their action on a basis vector |m) by 


(29) 


ipk\m) 

V'felm) 


lo 

jo 


-•} if mfc = o, 

otherwise, 

=•} if mfc = •, 

otherwise, 


where is the Maya diagram obtained from m by inverting the color of the 

marble on site k. Observe that the operators ifk and are adjoint to one another. 
In particular, i/fcV'fe (resp. V'fcV'fe) is the orthogonal projector on the space spanned 
by Maya diagrams m with mj, = • (resp. = o). Fermionic operators obey the 
following well-known canonical anticommutation relations: 


Proposition 9. For any k and k' in Z + 1/2, we have 

(30) = 0, = = 5k,k'- 

Flere {a, 6} denotes the anticommutator of a and b : {a, b} = ab + ba. 


Proof. Easy. 


□ 


Define the fermionic generating functions 

(31) -ipiz) = z^'ifk, = Y 

Proposition translates into 

(32) {il}{z),if{w)} = Q, {■!/*(z),i/>*(ui)} = 0, {'if{z),il)*(w)} = 5{z,w) 

where S(z,w) = J^kez+i(^/'’^)^ i® formal Dirac delta function. It is straight¬ 
forward to check that 

(0lif(z)if*(w)l0} = Y (^) = 1^1 > l^l> 

\w J z — w 

k^Z-\-lj2 ^ 

(33) 

{0\iI;*{w)'iIj{z)\0) = Y 

keZ+l/2 
k>0 

Here the leftmost equal signs correspond to formal identities, but the rightmost 
equal signs require to treat z and w as complex variables. Let us also mention a 
lesser-known fact about the action of the involution uj on the fermionic operators 
(recall that w is the involution that maps a charged partition (A,c) to (A/—c), 
hence can be seen as acting on F). 


\/ zw 
z — w 


for Iwl > Izl. 


Proposition 10. For k G Z + 1/2, we have 

(34) ujtfluj = (-I)<^+'=+i/V-fc, *-e. uj'if*{z)oj = 

where C is the charge operator (acting on the multiplication by c). 


Proof. Follows from the fact that, for any integer k' and any Maya diagram m of 
charge c, we have 

(35) #{j > k', mj = •} - #{j < fc', = o} = c - k'. □ 
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Last but not least, we have the following commutation relations between bosonic 
and fermionic operators. 


Proposition 11. Given two formal variables x,z we have 

(36) rL+ix)ip{z) = — ijj{z)TL+{x) 

1 — xz 

(37) rj{+(x)f’(z) = (1 + xz)f’(z)rji+(x) 

(38) rL+(x)-i/’*(z) = (1 - xz)-if*(z)rL+(x) 

(39) rji+(x)7f*(z) = —■if*(z)rji+(x) 

(40) i>{z)TL-(x) = (1 - x/z)Tl-(x)-iP{z) 

(41) fj{z)rii_{x) = —^—TR_{x)'tf(z) 

1 + XjZ 

(42) 'if*{z)TL_{x) = — ^-—rL_(x)tp*(z) 

1 — XjZ 

(43) tjj*{z)rii_(x) = (1 + x/z)rR_{x)ilj*{z) 

The first four (resp. last four) formal identities correspond to actual converging 
series when \z\ < x~^ (resp. \z\ > x). 

We give in Appendix a combinatorial proof of these identities (note that it is 
sufficient to establish only one of them, the others follow by taking duals, inverses 
and conjugates by ui). For algebraic derivations, see for instance the references 
given at the beginning of Section 


4.2. Constrained transfer matrices. The fermionic operators can be used to 
enumerate constrained dimer configurations. A hrst natural idea, already used 
in |OR03j . consists in inserting some orthogonal projectors fikifl. or i/'feV’fc (with 


various k’s) within the product of bosonic operators (22) forming the partition 
function, which has the effect of forcing black or white marbles to be present at 
given positions. However, this does not fully determine the positions of the dimers 
(there are ambiguities for the columns containing both horizontal and diagonal 
edges). Remarkably, for an arbitrary rail yard graph and an arbitrary finite set E 
of edges, there is a suitable way of inserting fermionic operators which precisely 
forces each edge of E to be covered by a dimer. 

We first introduce convenient notations. Recall that writing (a,/3) for an edge 
implies that its endpoints a and /3 are such that is even and is odd. Any 
finite set E of edges of a rail yard graph can be decomposed “column by column”, 
hence written in the form 

(44) E = i) , . . . , (cXi^rai 5 fii,mi)f (bl,l) ^i,l) ’ ’ ' ' ’ 5 )} 


where afj = 'yfj = 2z, = 2i — 1 and 6f j = 2z + 1. Here mi (resp. to') is the 

number of edges of E connecting vertices with abscissas 2z — 1 and 2i (resp. 2i and 
2z + 1), and is zero except for finitely many i. 

Theorem 12 (Constrained transfer matrix decomposition). Let E he an arbitrary 
finite subset of edges of the graph G = RYG(£, r, a, 6), which we decompose as in 


(44), and let Ui (i = i..r) be the number of diagonal edges of E with an endpoint of 
abscissa 2i. The sum of the weights 0 of all pure dimer configurations containing 
E reads 


( 45 ) 


W(G;az,A) = (0|r,T,+i.--r,|0) 
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where, for all i € [£--r], the constrained transfer matrix Ti is given by 


(46) = 


with ki = mi(mi — l)/2 + — l)/2 if Ui = R, and 




(47) = 


with ki = mi(mi — l)/2 + in^{m^ — l)/2 + Ui if at = L. More generally, the sum 
of the weights of all admissible dimer coverings with left boundary state 1 and right 
boundary state r containing E reads 

(48) W{G,\,v-x,E) = (l|r,T,+i..-r,|r). 


Observe that we recover (22) and (21) when E is empty. The order in which we 


take the products of fermionic operators in (46) and (47) is irrelevant, as long as 


we take the same order for both products from 1 to nii, and for both products from 
1 to m'i (otherwise, we might get a wrong sign). 


Proof of Theorem\lS\ It is sufficient to prove (48) in the case of an elementary RYG 


i.e. to prove that R is indeed the wanted constrained transfer matrix. The general 
case immediately follows by concatenation (i.e. by the transfer matrix method), as 
done in the proof of Proposition 

Let G = KYG{i,i,ai,bi) be an elementary RYG, 1, r two boundary states, and 
E a finite subset of edges of G. The general decomposition (44) reads here simply 


(49) E { (oi, /?x),..., (Oyy^ 5 jdm) J (T1 i ^ 1)5 ■ ■ ■ ; t ^rr 
Let US first assume that at = R, we then need to check that 

(50) IT(G,l,r;x„L;) = (l|T,|r). 


')}• 


with Ti given by (46). This is immediate in the case m' = 0: indeed the presence 


of a dimer on the edge {aj,j3j) (which is necessarily horizontal) is tantamount to 
having a white marble o at position aj = /3j in 1, see Figure 0A) . This can be 
achieved at the level of transfer matrices by multiplying the unconstrained transfer 
matrix Tai^bii^i) on the left by the projectors with k = 13^, j = l..m, and we 


get (46) upon anticommuting all to the left. 


The case m' ^ 0 is slightly more involved and requires the introduction of suitable 
“particle hopping” operators. Recall that the notation denotes the Maya 
diagram obtained from m by inverting the color of the marble on site k. Let us 
consider an edge Cj = (pjj,6j)'. having this edge covered by a dimer implies that 
Ly = r^y = •, but the converse is not true. However there is a bijection between, on 
the one hand, admissible dimer configurations with boundary states 1 , r containing 
Cj and, on the other hand, admissible dimer configurations with boundary states 
l( 7 j that necessarily contain the edge e' on the left of 7 ^, see Figure 9lB- 

C). We deduce that, at the transfer matrix level, we can force the presence ot Cj 
by multiplying the unconstrained transfer matrix Fq. b (xi) by on the left and 
by ipgy on the right (if 1 or r do not have black marbles at the required positions 

^3 

then (1| or |r) will be killed by this ip or ip*, as it should be). More generally, the 
product (l|V’ 7 f ■ ■ ■ V’ 7 ^ ai,bi(xi)il)*gy ■ ■ ■ 'pgy\'r) is nonzero if and only if there is a 
dimer configuration with boundary states 1, r containing all the ej, j = In 

that case this configuration is unique and the product is equal to xf, with n the 
number of dimers on diagonal edges other than the Cj (it is easily seen that the signs 
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Figure 9. Localization of dimers in the case of the elementary 
graph or type R— (the case of i?+ is similar). Marbles (resp. 
edges) whose status is not fixed by the discussion are represented 
with question marks (resp. dotted lines). (A) There is a horizon¬ 
tal dimer in the left column if and only if the left odd vertex is 
occupied by a white marble. This marble is localized by applying 
the operator (B) Configurations with a horizontal dimer in 

the right column are such that both marbles on this level are black 
(although this condition is not sufficient). They are in bijection 
with configurations with a dimer in the left column such that both 
marbles on this level are white. The operator (resp. 'ipk) in¬ 
serted on the right (resp. left) has the double effect of switching 
the color of marbles on the fc-th level and of forcing the colors of 
these marbles. (C) The case of a diagonal dimer. 


induced by the '0/V'* all cancel out). Upon reordering the '0*’s in reverse order, 
multiplying by a;"* (to account for the weight of the dimers on diagonal edges in 
Cj, j = l..m') and multiplying by projectors on the left, we conclude that (50) is 
true with Ti given by (461. 

The case ai = L can be deduced from the discussion of the case Qi = R, by 
performing a central symmetry and exchanging the colors of the marbles. In other 
words, we simply need to take the dual of (46) (vertical symmetry) and conjugate 
with the involution uj acting on charged partitions/Maya diagrams (horizontal and 
color symmetry). There is a slight subtlety regarding the sign though, which can 
be treated using Proposition When taking the dual of (46), the order of the 
operators is reversed, rfl;_±(a;i) is changed into and each ip is changed into 

a Ip* (and vice versa). When conjugating by w, ri^^zp(a;i) is changed into rL_zp(a;i) 
and each ip* is changed back into a ip (and vice versa), up to a sign. By ( [M| , the 
signs cancel out for horizontal edges in A, but combine into a sign —1 for each 
diagonal edge in E, which explains why the sign (—1)^’ is different in (47). □ 
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Remark 13. Similar bijective arguments are used in Appendix to prove the 
commutation relations between bosonic and fermionic operators stated in Proposi¬ 
tion [TTJ 


4.3. Prom the Schrddinger to the Heisenberg picture. Theorem 12 expresses 
W (G; X, E) as a product of bosonic and fermionic operators taken between two 
vacuum states, which we may rewrite using a strategy coming from |OR03j . similar 
to that used in Section 3.3 for the proof of Theorem]^ move all P+’s to the right, 
and all P-’s to the left, so that they are absorbed by the vacuum states at the end. 
In this process, we first pick multiplicative factors due to the commutations between 
P+’s and r_’s: those are precisely the same as for the partition function Z{G',x). 
Second, the fermionic operators get “conjugated” by the P’s crossing them. All this 
allows to rewrite 

(51) W{G;x,E) = Z{G-,x){mTe+i ■ • •T;|0) 

where, for all i G [£..r], we set 

(52) 

( mi \ / rui \ / \ / 

with, for /?, a respectively odd and even vertices of G, 

( \ 




n ra.,+ (a;i) Yl 


(53) 


V 


z>rr/2i 

bi = - 




/ 


'l'(a) = Ad 


n 


- 


{x^) 


n 


ra.!-(^.) 


z<q;^/ 2 if ai—L 
j2 if ai — R 


i'>a^/2 if ai—R 
2>a^/2 if ai—L 
bi^- 


ll^a 




Here Ad denotes the adjoint action: 

(54) Ad{A)B = ABA-^ 

with A, B operators acting on E, with A invertible (recall that, by ( [l7| , this is the 
case for the P’s). Two remarks are in order. First, we may put the P+’s and the 
P_’s in any order we want in (53), as this does not change their adjoint action. 


Second, by Proposition 11 'I'*(/3) and 'l'(a) are (formal) linear combinations of 
V'*’s and '0’s, respectively. 


Remark 14. In physical terms, passing from the 0*/0’s to the 4'*/4>’s can indeed 
be interpreted as going from the Schrddinger to the Heisenberg picture of quantum 
mechanics (the abscissa playing the role of time). It appears from (53) that creation 
and annihilation operators are naturally attached to respectively even and odd 
sites. An analogous situation appears in the path integral formalism: if K denotes 
the Kasteleyn matrix of a finite planar bipartite graph, then the determinant of 
K, yielding the dimer partition function, can be written as a Grassmann-Berezin 
integral |ID89) 


(55) det a: = J 

i 3 

where the ^i’s and are Grassmann variables attached to the even and odd vertices 
of the graph, respectively. Furthermore, the contribution of dimer configurations 
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containing a given collection of edges (zi, ji),..., {is,js) is proportional to 


(56) 



•aCj e 






In other words, dimer correlations are given by the expectation value of a product 
of fermionic operators, whose form is reminiscent of (52). Note however that the 
approach followed in this paper is more akin to canonical quantization. 


4.4. 


Proof of Theorem By (51) and (52), the ratio 


(57) 


Pg-,AE) 


W{G;x,E) 

Z{G]x) 


is given by a product of fermionic operators taken between two vacuum states. This 
can be rewritten as a determinant using Wick’s formula, as follows. 

For an even vertex a and an odd vertex /3, define the naturally ordered product 
(or “time-ordered product”) of 4'(a) and 'I'*(,5) by 


(58) 


r(vI/(a),vI'*(/3)) 


4'(a)«'*(/3) 


We may more generally consider the naturally ordered product of more than two iP 
or dt*, by ordering them according to the abscissa of their argument and multiplying 
by the sign of the corresponding permutation (note that operators with the same 
abscissa anticommute, hence their order is irrelevant). Observe that ■■ - Tr 

is, up to a factor, the naturally ordered product of the if and associated with 
the endpoints of the edges of E. Denoting now by {ei,..., Cg} the edges of E, 
with Ci = and by H{E) the number of horizontal edges in E whose right 

endpoint is at an even abscissa, we have 


Pa-dE) 

(59) 





(0|r ('k(ai), 4/*(A),..., 4'(a«), 4'*(/3,)) |0) 




det (0|r(4-(aO,4/*(^,))|0) 


where we apply Wick’s formula to pass from the first to the second line, and where 
the sign (—1)-^*^^) arises from the reordering of the fermionic operators (in partic¬ 
ular, dimers having their right endpoint at an even abscissa appear in the “wrong 
order” in the naturally ordered product, but the resulting sign is cancelled in the 
case of diagonal dimers by that present in Theorem |12[ ). For completeness, we pro¬ 
vide a detailed derivation of Wick’s formula in Appendix |B| To complete the proof 
of Theorem it remains to check that 

(60) C„.^ = (0|r(4'(a),4/*(/3))|0) 

has the announced expression ( [To| ) . At this stage we need to discuss a bit analyticity 
conditions (so far all our computations were done by treating the weights x^’s 
as formal variables). Recall that we aim at proving Theorem under the mere 
assumption that the partition function Z{G;x) is a convergent sum, which boils 
down to the conditions ([|) and ([rWsee again Remark [§. 

Let us temporarily strengthen (M into the condition 

XiXj < 1 for all i < j such that bi = + and bj = — 


(61) 
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(that is, we also impose XiXj < 1 when ^ aj). We introduce the quantities 

1 


(62) 


p+{k) = inf 
2i<.k 


P-(k) = sup Xi 
2j>k 


0 ’ 


which are nonnegative nonincreasing functions of k, such that P-{k) < p+{k) by 
(61) and ([^. Recalling the definition © of the fermionic generating functions 
'tp{z) and ip*{z), the definition (53) of 'l'(a) and 'I'*(/3), and the bosonic-fermionic 
commutation relations of Proposition El we may write 


(63) 




Fci^{z)'ip{z) dz 




<« = 


dw 


/C^ Z” z Z'ivi Fpx{w) w 

with Cz (resp. Cyj) a circle centered at 0 with radius comprised between p_{a^) and 
p+{cF) (resp. between p-{f3^) and p+{(3^)). We deduce that (0|T ('I'(a!), 'l'*(/3)) |0) 
is equal to the wanted expression (lIo|) provided that 


(64) 


(0|r {ip{z),tp*{w)) |0) = 


/zw 


But this readily follows from the definition ( |58[ ) of the naturally ordered product 
and from the relations (33), provided that we take the radius of Cz to be strictly 
larger than that of Cw if and vice versa otherwise. We may now freely 

deform the contours Cz and C^,, as long as we hit no pole of the integrand; this 


establishes (10) under the conditions (i)-(iii) for the contours, hence Theorem]^ 
under the assumption ( |61[ ). 

We now explain how to relax this assumption into (§. We proceed by multiplying 
each Xi by a factor t G [0,1], and noting that the assumption (61) hence the identity 
(11) are satisfied for t small enough. We will show that both sides of (11) have an 
analytic continuation in < to a domain containing the closed unit disk. As apparent 
from Theorem this is the case for the partition function Z{G',tx) (provided that 
^ is satisfied, of course). The quantity W{G\ tx, E), being a sum over a restricted 
subset of configurations, is analytic too, hence so does PG-tx{E) which is the left- 
hand side of (11). As for the right-hand side, let us show that we may find two 


contours Cz and such that as defined by (10) with the XiS multiplied by 


t, is manifestly analytic as t varies over the unit disk. For t = 1, let us introduce 
the quantities 


(65) 


PR+{k) = inf 

2i<.k 


PL+{k) = mf 

*:(ai,bi)=(L,+) 

2i<k 


Xi 


PL-{k) = 


rP 


PR-{k) = 


j:{a. 


sup Xj 

,bj)={L-) 
2j>k 

sup Xj 

M={R-) 

2j>k 


which are refinements of p+{k) and p-{k). We take Cz to be a circle centered 
on the real axis, surrounding the real interval [—pij_(a’'), 0] but not intersecting 
[pL+icF),+oo[. Similarly, we take Cw to be a circle centered on the real axis, 
surrounding the real interval [Q, pL-{j3^)] but not intersecting ] — oo,—pi^+(/?’')]. 
We furthermore want Cz and C^ not to intersect, with Cz surrounding C^ iff 
this is possible because the functions pl± are decreasing with k and 
such that pL-{k) < pL+{k) for all k hence, for < /3^, we have pL-{P^) < 
PL-{a^) < pl+{o:^) so that the “window” ]pi_(/l^),pL+(a^)[ in which both Cz 
and Cw must pass is nonempty, see Figure 10 (and the case a® > is treated 


by a similar reasoning on pr±). Observe that these contours satisfy precisely the 


conditions (i)-(iii) stated below equation (10). Keeping these contours fixed, we 
let t vary in the interval [0,1]: no pole of p^{z,w) ever hits the contours, the 
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Figure 10. The contours Cz and Cu, in the case > a^. The 
upper (resp. lower) dashed regions are the intervals that the con¬ 
tour Cz (resp. Cw) must avoid. On this figure it is assumed that 
PR-{a^) > but the converse can also hold. 


conditions (i)-(iii) remain satisfied, and we conclude that (10) defines the wanted 


analytic continuation of 


□ 


4.5. Inverse Kasteleyn matrix. Let us now connect Theorem]^ to the general 
Kasteleyn theory for the dimer model on planar bipartite graphs. We mention that 
a similar connection was already observed in |OR07j for the case of lozenge tilings 
corresponding to skew plane partitions, see also [BF141 Section 5]. 


Definition 15. Let G = (V, E) be a planar bipartite graph with no multiple edges. 
Consider a collection of weights on the edges uj '■ E ^ IR_|_. 

A Kasteleyn orientation of G is a map p : E ^ 1} such that for any face E 

the product over the edges surrounding E gives : 


( 66 ) 


n 

eeF 


1 if F is of degree 2 (mod 4) 
— 1 if F is of degree 0 (mod 4) 


The Kasteleyn matrix of G associated to p is the matrix K, whose rows (resp. 
columns) are indexed by white (resp. black) vertices of G, such that for any couple 
{w, b) of a white vertex and a black vertex. 


(67) 


JC{w, b) 


0 if w not adjacent to b 

p{w, b)uj{w, b) if there is an edge {w, b) 


Local statistics for dimers are known to be given by determinants of submatrices 
of the Kasteleyn matrix |Kas61[ IPer691 IKen97| : 


Theorem 16. Let{G,uj) be a finite weighted planar bipartite graph. The probability 
that the dimers ei = (ici, 6i),..., = (ws^bs) are present in a random dimer 

configuration sampled with a probability proportional to its weight is: 


( 68 ) 


P{ei,...,es) = j A(wi,&i) j det ^{b„Wj)^ 




Let us go back to the rail yard graphs. In this case, the black (resp. white) 
vertices are the even (resp. odd) vertices. Recall that in pure dimer coverings of 
a rail yard graph G , all the odd vertices with negative (resp. positive) ordinate 
on the left (resp. right) boundary are unmatched. These pure dimer coverings on 
G correspond to classical dimer coverings with finitely many diagonal edges on 
the graph G, where those unmatched vertices and the edges attached to them are 
removed. Compare Figure 11 and Figure]^ 
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Figure 11. The perfect matching of the modified rail graph cor¬ 
responding to the pure dimer covering of Figure]^ 


Let us denote by 9JI the set of matched odd vertices of G. Define rj E ^ {~li 1} 
by 

1—1 if e is a horizontal edge whose right end is an even vertex, 
(69) 77 (e) = <^ . 

1 otherwise. 


It is easy to check that 77 is a Kasteleyn orientation on G (here the faces have degree 
4, 6 or 8). Construct the infinite Kasteleyn matrix JC on that graph as in (671. The 
restriction /C of /C to rows indexed by 911 is a Kasteleyn matrix for G. We now 
relate our correlation kernel C to the infinite Kasteleyn matrix /C: 


Theorem 17. Let G = RYG(£, r, a, &) be a rail yard graph and 1C be its Kasteleyn 
matrix for the previously defined orientation. Then if C is the matrix defined in 
Theorem^ we have: 

(1) for any a, a' even vertices, {CIC)a,a> = ha,a',' 

(2) for any P,/3' odd vertices in DJI, {lCC)p,j 3 i = d/ 3 ,/ 3 '. 


We state a lemma that will be useful to prove the theorem. Recall that we 
have a bosonic operator Tji±{xm) (resp. TL±{xm)) at position (2 to-|- 1/2,0) (resp. 
{2m — 1/2,0)) when = R (resp. = L)- To simplify notations, we will also 
place a bosonic operator Id at every position {2m — 1/2,0) (resp. {2m + 1/2,0)) 
when Um = R (resp. Om = L). This does not change the naturally ordered product 
of operators, and now we have one bosonic operator at each half-integer abscissa in 
[2£ — 1,..., 2r + 1]. Let B~ (resp. B^) be the bosonic operator at abscissa i — \ 
(resp. i + |). Denote by x ~ 77 the fact that two vertices x and y are adjacent in 
G. 

Lemma 18. We have the following properties: 
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(1) Let a be an even vertex at position {i,k). Then 

Y, IC{P,a)rpyB- = -B-rk, 


(70) 

(71) ^/C(/3,a)B+V'^. =i/.^i3+. 
(2) Let P be an odd vertex at position {i,k). Then 


^r^OL 


(72) 

(73) 


^/C(/3,a)V'ay5 =B ijjk] 


CX.r^^ 

0'^<i 


Y IC{j3,a)B+tjjay = -ipkB^. 


a^(3 

a^>i 


Proof. We will prove only identity (70). The other identities can be proved in a 
similar way. We distinguish three cases : 

Case 1: B~ = Id. To its left, a has only one neighbour /3, which is at height k. 
We conclude from the fact that /C(/3, a) = —1. 

Case 2 : B~ = Tl-^x). To its left, a has one neighbour fdi at height k and one 
neighbour /?2 at height k — 1. We compute 

Y )^iP,aW^y^L-{x) = JC{/3i,a)iplTL-ix) +/C(/32,a)'0LirL-(a^) 


(3^CK 

'^^<i 




^-k 


-k 


^-k 


^-k 


if*{z)rL-{x) + 
fj*{z)rL-{x) + 


,-fc+i 


xif*{z)TL-{x) 


,-k 


-r{z)TL.ix) 


i--]r{z)rL-{x) 


rL-{xW{z) 


= -rL_(a;)i/>, 


k 5 


where we used proposition 11 to switch 'ip*{z) and T l-{x). 

Case 3 : B~ = T l+{x). In the proof of case 2, replace T l-{x) by rL+(a:), k — 1 
by fc + 1 and | by xz. □ 


We can now prove Theorem 17 


Proof of Theorem lPll In this proof, to make notations lighter, Z{G]x) will be ab¬ 
breviated as Z. Let us prove the first part of the theorem. Fix two even vertices a 
at position {i,k) and a' at position [i',k'). Then : 


(74) 


(C^)(a',a) = ^C(a',/3)/C(/3,a). 


This sum has at most three nonzero terms, corresponding to the three odd neigh¬ 
bours of a, at the abscissas i — 1 and i -|- 1. We now distinguish according to the 
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position of i' relatively to i. 

Case 1 : i' < i. 

If /3 ~ q; is at abscissa i — 1, C{a', /3) can be written in the following form : 

(75) c(a',/?) = 

where and F^^^ are the products of the bosonic operators located respec¬ 

tively before the abscissa i', between the abscissas i' and i — 1 and after the abscissa 
i -f 1. 

Similarly, if ~ a is at abscissa i + 1, C(a',P) can be written in the following 
form: 

(76) C(a',/3) = |(0|F(i)^fc,F(2)s-B+7)>^,F(3)|0). 

For each j3 a appearing in the sum of equation ( |74| ) , we are going to move the 
'tp*py between B~ and . We separate the cases j3^ <i and (3^ > iio apply part 1 
of the lemma: 

(C/C)(a',a) = ^/C(/3,a)C(a',/3) 

/3~a 

= ^ /C(/3,a)C(a',/3) -f ^ JC{l3,a)C(a\ P) 

jSr^a ^'^Oi 

Z/ z 

= 0 . 

Case 2 : i' > i. 

Here ipk' is to the left of any so each C(a',/3) comes with a minus sign. 
Using here again part 1 of the lemma to move tp^y between B~ and B^ for all 
three terms, we get similarly that (CJC){a', a) = 0 in this case. 

Case 3 : i' = i. 

If (3 is at abscissa i — 1, then C(a', /?) can be written in the following form: 

(77) C(a',/3) = -^(0|F(i)V'^,H-V'fc,i?+F(2)|0), 

where the minus sign comes from the fact that ip*py appears before i/'fe'- 
If /3 ^ a is at abscissa i -f 1, we have: 

(78) C{a',l3) = ^(0|F(i)H-^fc,H+^^,F(2)|0). 

Applying again part I of the lemma, we get: 

{CJC)(a', a) = Y, ^(/3. a)C(a\ /3) 

/3~a 

= ^ JC{(3,a)C{a\(3) + ^ IC{P,a)C(a', (3) 

= 4(0|r(^)r(2)H-^/:^V'fc'-B+r(3)|0) + l(0|F(i)F(2)H-^fe,^*H+F(3)|0). 

z z 

Using that tpl'ipk' + V'fe'V'fe = (Proposition]^ and that 

(79) Z = (0|f(i)f(2)b-H+f(3)|0), 
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We conclude that in this third case, {CIC){a',a) = 6k,k'- This concludes the 
proof of part 1 of the theorem : we have shown that for any even vertices a and 
a', (C/C)(a',a) = 6a\a- 

To prove part 2 of the theorem, note that if f3 and /?' are two odd vertices in 9Jl, 
then 

(80) (/CC)(/3, /?') = ^ /C(/3, a)C(a, /3'). 

a~/3 

We will only treat the case when (3 is on the left boundary. The case when (3 
is on the right boundary is similar to this one. The case when (3 is in the bulk is 
similar to the proof of part 1 of the theorem, by making use this time of part 2 of 
the lemma. Let (2£ — 1, k) be the coordinates of /3 and (i', k') be the coordinates 
of (3'. Since (3 G SEJl, we have k > 0. Note that /3 has neighbours only to its right. 
Again we distinguish according to the position of f3'. 

Case 1 : i' >2i- 1. 

If a ~ /3, C(a, (3') can be written in the following form : 

(81) C(a,/30 = ^(0|i3+i/>„.rT)V;^r(2)|0). 

Using the second equation of part 2 of the lemma, we obtain: 

(/CC)(/3,/3') = 5]/C(/3,a)C(a,^') 

q~/3 

= ^ JC{(3,a)C{a,l3') 

a~/3 

a^>2£-l 

= -^(0|V>fci?+r«V;*,r(2)|0). 

Since A: > 0, we have {0\tpk = 0. So in this case, (/CC)(/3,/?') = 0. 

Case 2 : i' = 2i- 1. 

If Of ~ /3, C(a, /3') can be written in the following form : 

(82) C{a,(3') = -l(0|V:^,i3+V:„.r(i)|0). 

Using the second equation of part 2 of the lemma, we obtain : 

(83) (/CC)(/3,/3') = +i(0|V-:,V’fcB+r(i)|0). 

Using Proposition]^ we get : 

(/cc)(/?, /?') = -1 m^krk'B+r^"^ |0) + 6k,k' I (0|i?+r(i) |0) 

= 0 + 6k,k'~^Z = 6k,k' ■ 

To sum up, if /3 is an odd vertex in 3Jl on the left boundary and j3' is any other 
odd vertex in 9H, (lCC)i^p,pi) = 6p,pi. □ 

Remark 19. If C (resp. K.) is the restriction of C (resp. K.) to columns (resp. rows) 
indexed by the set 9J1 of matched odd vertices, then C is a left and right inverse of 
JC. Indeed, the extra terms coming from vertices not in DJI vanish when evaluated 
against left and right vacuums. 


28 


C. BOUTILLIER, J. BOUTTIER, G. CHAPUY, S. CORTEEL, AND S. RAMASSAMY 


5. Lozenge and domino tilings 


In this section we discuss how the RYG dimer model gives rise, upon taking a 
constant or an alternating Li?-sequence and performing a simple change of coordi¬ 
nates, to a class of lozenge or domino tilings of the plane, respectively. The class 
of lozenge tilings, discussed in Section |5.1[ contains in particular plane partitions 
and was already studied in |OR031 IOR,07j . The class of domino tilings, discussed 
in Section 5.2 is that of steep tilings [BCC14| . and contains tilings of the Aztec 
diamond and pyramid partitions as special cases. 

5.1. Lozenge tilings. In this section we consider the case where the Li?-sequence a 
is constant. We assume that a = (it is easy to check that taking a = 

gives rise to the same model up to a vertical reflection and an inversion of the sign 
sequence). In this case the only two elementary RYG involved are of type L+ 
and L—. According to the discussion of Section the transfer matrix operators 
corresponding to these two graphs can be interpreted as the operators ri_|_ and 
r, whose action on the Bosonic Fock space interlaces a partition “upwards” 


or “downwards” respectively, see (14). It follows that, for each sign sequence 
b G {-f, —admissible dimer coverings of RYG(£,r, a,6) are in bijection with 
sequences (A(*))^<i<r+i of integer partitions such that where the inter¬ 

lacing relation is ^ or if 6^ = -|- or 6^ = —, respectively. The correspondence goes 
via Maya diagrams and is the one described in Section It is well-known that such 
sequences are in bijection with certain lozenge tilings of the plane, see e.g. [OR03| . 
so the reader may be in familiar ground. In the rest of Section |5.H we will just 
sketch how to recover those lozenge tilings from the rail yard graphs. 

We first apply the following coordinate transformation to each vertex of the rail 
yard graph RYG(£, r, o,6): 


/ 


(84) (2£-l+x,g) 


2£-l + -[x/2] 


V 


V5, ,01 

— L^/2j,y-- ^ 
2=1 


(-l)L 


Figure displays the effect of this transformation on the elementary graphs of 
type L+ and L—. The transformation is designed in such a way that after the 
transformation all angles between incident edges are equal to ^, and that all edges 

have equal length (equal to ^). Therefore the concatenation of these graphs 
generate a portion of the regular hexagonal lattice, see Figure [Mi-Left. Since the 
planar dual of the hexagonal lattice is a triangular lattice, any dimer covering of 
the hexagonal lattice induces a covering of this triangular lattice by lozenges (each 
dimer connects two vertices in the primal, that correspond to two triangles in the 
dual, and the union of these two triangles forms a lozenge). See Figure [l4]-Right . 
We thus recover the promised class of lozenge tilings of the plane. 

Note that the fundamental covering of RYG(£, r, a, b) projects (via the coordinate 
transformation (84) and the dualization to lozenges) to a lozenge tiling in which 
all lozenges under a certain separating path are of “horizontal” type, while all 
lozenges above this path are of one of the two “vertical” types, see Figure This 
separating path is the image under the transformation (84) of the horizontal axis in 


the original embedding of RYG(£, r, a, 6). Equivalently, this path is a lattice path 
taking up (-I-) or down (-) steps, encoded by the sequence b. This path and its 
image are represented by red dotted paths on the figures of this section. 

Finally, note that Figure can naturally be seen as a 3-dimensional picture, 
namely as a portion of the boundary of the region (K^\A) xK+, where A is the shape 
of an integer partition. Here the contour of the shape A coincides with the separating 
path just mentioned, i.e. the partition A is encoded in Russian notation by the 
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O ^ 

■ -^+ ■ -^+ ■ L_ '■ L+ : L_ '■ L_ 


Figure 12. A rail yard graph with a constant LR sequence a = L° 


and sign sequence b = + H-1— 

■ O 


Q.■ O Q 



>o 

O.■ o o 


L+ 

Figure 13. Deformations of the 
L— that generate a “honeycomb’ 


, equipped with a pure covering. 


O .o o 

• o 



elementary RYG of type L+ and 
’ lattice by concatenation. 


sequence b, see Figure again. Note also that all the lozenge tilings corresponding 
to pure dimer coverings of RYG(£, r, a, ^ are obtained by “adding cubes” to this 
3-dimensional diagram in such a way that the heights of cubes stay nonincreasing 
on horizontal coordinates, going away from the axes. In particular the case b = 
_l_n_n corresponds to plane partitions of half-width at most n. This 3-dimensional 
interpretation is well known and is the one already considered in [OR 03) . 
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Figure 14. Left: The dimer covering of Figure 12 as a covering 
of the honeycomb lattice, via the transformations of Figure 
Right: The same objects displayed as a “lozenge tiling” of the 
plane. 



Figure 15. The fundamental lozenge tiling corresponding to the 
case a = L^ and 6 = + H-1-. 
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5.2. Domino tilings. We now consider the case where the sequence a has even 
length and is an alternation of L and R. Up to an elementary symmetry we can 
assume that a = [LR)^, where r —£+1 = 2k, and we fix an arbitrary sign sequence 
b G It follows from a simple inspection of the face types of elementary 

rail yard graphs (Figure |^a)) that all the inner faces of the graph RYG(£,r,a,b) 
have degree 4 or 8. Moreover, for 0 < f < fc, odd vertices located at the abscissa 
X = —2£ — 1 + 4i + 2 all have degree 2 (since they lie at the interface between an 
L-type and i?-type elementary graphs, from left to right). These degree 2 vertices 
are bounded by faces of degree 8 on their two sides. See Figure [T6| for an example. 

We now let R(i, r, a, b) be the graph obtained by contracting all the edges incident 
to these inner vertices of degree 2: in this graph all inner vertices of degree 2 have 
disappeared, and all the non-boundary faces have degree 4, see Figure [T7|-Left for 
an example. Moreover, it is easy to see that each dimer covering of RYG(£, r, a, 6) 
induces a dimer covering of R(£,r,a,b), with the same boundary conditions (just 
forget dimers on contracted edges and leave the other dimers as they were). The 
coordinate transformation on non deleted vertices that goes from RYG(£, r, o, 6) to 
R{£, r, a, b) is given by: 


(85) 


(-2£- l + x,y) 


-2i-l + x-2 


X 2 


,y 


Since all the inner vertices and faces in this new graph have degree 4, this new 
graph is isomorphic to a portion of the square lattice. We let the reader check that 


this isomorphism can be made explicit by composing (85) with the transformation: 

( 86 ) 


(-2£ -l + x,y) I—)■ {-2e -1+x + y- K,,, y - K ,,), 


where Figure 17 for an explicit example. The image of 


R{(., r, a, b) via these transformations is a portion of the square lattice lying in the 
oblique strip {(X, Y) : |Y — X -|- 2£ -|-1| < 2k}, see again Figure 17 

Similarly as in Section |5.1[ we can dualize this picture to switch between a 
description in terms of dimers to one in terms of tilings. The dual of the square 
lattice is again the square lattice, and any dimer in the primal induces a domino (the 
union of two adjacent squares) in the dual. We thus recover a model of tilings of the 
oblique strip by dominos, which are exactly the steep tilings introduced in [BGG14| . 
We invite the reader to consult this reference for a thorough discussion on steep 
tilings, their link with height functions, their encoding in terms of partitions, etc. 
Here we just mention again that, not only do we recover here the enumerative 
results already proved in |BCC14] , but we obtain the inverse Kasteleyn matrix and 
dimer correlations for these models, up to the changes of coordinates described 
above. 

Two particular subclasses of steep tilings had been considered previously. The 
first one is the class of domino tilings of the Aztec diamond, which corresponds to 
the case where the sequence b is also alternating. As far as we know, this is the 
only case for which the inverse Kasteleyn matrix has been computed before |Hel001 
IGY14] . The other one is given by pyramid partitions |You09] . that correspond to 
the case where b = See [BGG141 Section 4.2] for their connection with 

steep tilings. We leave as an exercise the task of making fully explicit the changes 
of coordinates and the calculations for pyramid partitions, similarly as we will do 
for the Aztec diamond in the next section. 


Remark 20. Going from rail yard graphs to steep tilings induces, in a sense, no 
loss of generality. Indeed, by forbidding diagonal dimers in the column z of a rail 
yard graph (by taking their weight Xi to be zero), this column becomes trivial 
in the sense that its two boundary states are necessarily equal, and hence it can 
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Figure 16. A rail yard graph with LR sequence a = {LR)^ and 
sign sequence b = (H—)" for n = 4, equipped with its fundamental 
dimer covering. 


be contracted in the graph (i.e. we can drop the corresponding elements from the 
LR and sign sequences). Starting from an infinite alternating LR sequence, it is 
possible to produce any LR sequence by performing such contractions, and hence 
to obtain any RYG. The coordinates of rail yard graphs however allow to express 
the dimer correlations in a more compact form, which is one of our motivations for 
introducing them. 


6. The Aztec diamond revisited 


In order to illustrate Theorem on a concrete example, we concentrate in this 
section on the particular case of domino tilings of the Aztec diamond of size 
n |EKLP92al lEKLPQ^ , which can be obtained as a rail yard graph associated 
with the sequences a = (LR)^ and b = (H—)", see Figures 16 and 17 We will 
suppose that £ = 0 and r = 2n — 1, and denote by the corresponding rail yard 
graph. Note that, though Gn is an infinite graph, its admits only a finite number 
(2"("+i)/^) of pure dimer coverings. Those pure dimer coverings coincide with the 
fundamental covering outside a finite region which is isomorphic to the Aztec dia¬ 
mond graph of size n, see |BCC14l Section 4.1] for a discussion of this phenomenon 
in the language of steep tilings. The Aztec diamonds are essentially the only RYGs 
having a finite number of pure dimer configurations (this can be seen from the 
enumerative results of Section 2.4 and Remark 20). 

Let us now discuss the probability distributions over domino tilings of the Aztec 
diamond that we are considering. Recall the definition Q of the weight of a config¬ 
uration in the multivariate RYG dimer model. For a generic sequence xq, ..., X 2 n-i, 
we obtain the so-called Stanley weighting scheme |Pro971 IYan91| , see also |BGG141 
Remark 2]. The partition function reads 


(87) 


^(Gn,x)= (1+XiXj) 

0<i<j<2n—l 
i even, j odd 


which is a polynomial in the a:’s. By specialization we obtain the following distri¬ 
butions considered originally in [EKLP92a] : 

• the uniform distribution, obtained by taking Xi = 1 for all i. 
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Figure 17. Left: The graph obtained by contracting all the ver¬ 
tices of degree 2 in the RYG of Figure [Tbl that is, all the vertices 
of abscissa congruent to 2 mod 4. Right: the image of the graph 
on the left by the linear transformation i—)■ The graph 

is a portion of a square lattice of mesh 1. On both sides, the 
region of the graph where a covering may differ from the funda¬ 
mental one is represented in red. One recognizes the shape of the 
Aztec diamond of size n = 4. Note that the coordinate transfor¬ 
mation on uncontracted vertices to go directly from the graph of 
Figure 16 to the graph on the right is given by 


where (/)(x) = x — 2[^^J if a: > 0 and <^(—a:) = —(j){x). 


• the biased distribution, obtained by taking Xi = 1 for i even and = A > 0 
for i odd (or equivalently Xi = a/A for all i): this corresponds to attaching 
a weight A to each pair of diagonal dimers (which become vertical dominos 
in the Aztec diamond picture), 

• the so-called g™' distribution, obtained by taking Xi = g* for i odd and 
Xi = g“* for i even (which is the g-RYG specialization), and more generally 
the biased g''°^ distribution, obtained by taking Xi = Ag® for i odd and 
Xi = g“® for i even. 

Let us mention that our present approach does not seem to apply to the two- 
periodic weighting considered in |GY14[ Section 6], nor to the weightings considered 
in |DFSG14| . 

The study of correlations in domino tilings of the Aztec diamond has been a 
popular topic (especially among the members of the “domino forum”) and there are 
many previously known results, published or unpublished. In the rest of this section, 
we rederive several of these results as consequences of our general Theorem 


6.1. The biased creation rate and edge-probability generating function. 

Let us consider the Aztec diamond of size n in the natural coordinates with the 
origin at the center. We are interested in the probability to find a domino of a 
given type at a given position, under the biased distribution. Recall that we may 
distinguish four types of dominos: north-, south-, west- and east-going [GEP96| . 
By symmetry it is enough to consider only one type of domino, and we denote by 
P\{x, y, n) the probability that {x — 1/2, y) is the center of a west-going domino in 
a biased random tiling of the Aztec diamond of size n. Due to parity constraints, 
this probability vanishes unless x and y are integers such that a: -I- y -I- n is odd. 
As apparent from Figures lb and[^ west-going dominos correspond in the RYG 
setting to diagonal dimers in columns of type L+ (which makes them easier than 
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north-going dominos to deal with, since there can be no spurrious diagonal dimers 
outside the “interesting” region). 

Getting an expression for Px(x,y,n) (or its analogues for other domino orienta¬ 
tions) ameneable to asymptotic analysis has been of interest to several people, for 
instance it is used in [CEP96| as a way to proving the arctic circle theorem and its 
generalization to arbitrary A. However, this paper states the required expressions 
without proofs, and refers instead to a preprint by Gessel, lonescu and Propp that 
has not appeared so far. The proof for A = 1 (uniform distribution) can be found 
in Helfgott’s senior thesis |HelOO| . More recently, an expression for the related gen¬ 
erating function (still for A = 1 only) was proved by two methods in Du’s master 
thesis [Dullj . and was used in [BPllj as yet another route to the arctic circle the¬ 
orem (we note that those two references mention another lost “DGIP” preprint). 
At the suggestion of James Propp, which we thank for pointing out this gap in 
the literature, we will explain how expressions for P\{x, y, n) (for general A) can be 
obtained as applications of our Theorem 

The first step consists in going from dominos to dimers on RYGs. Using the iden¬ 
tification discussed above, a west-going domino centered on (cc —1/2, y) corresponds 
to a dimer on the edge (a, /3) with 

(88) a = (2m,y-l/2), /3 = (2m-l,y-kl/2), m = g [l..n]. 

By Theorem we immediately deduce the expression 


(89) 


Px{x,y,n) = 



m 


r i + \/w y+^-'^ 

\ 1 -I- A/z 7 




{w/z)y 


—-dzdw 
w) 


where Cz is a positively oriented contour containing 0 and —A in its interior, but 
not 1, and is a positively oriented contour containing Cz in its interior. Simpler 
expressions can be obtained for two related quantities: the so-called biased creation 
rate 

(90) CrA(a;, y, n) = {P\{x, y, n) - Px{x Pl,y,n- 1)) 

and the edge-probability generating function 


(91) 


Ilxiu,v,t) = } Px{x,y,n)u^vyc 


x,y,n 


Note that, in [GEP96| . the bias p is related to our A by p = 1/(1 -I- A), and the 
biased creation rate is expressed in terms of north-going dominos hence the present 
definition is adapted to the case of west-going dominos. 

Let us first consider the biased creation rate. Taking the difference Px{x, y, n) — 
Px(x -I-1, y, n — 1) in the double contour integral (891, the integrand is multiplied by 
a factor 1 — leading to a cancellation of the denominator {z — w), 

hence to the factorization 


(92) 


CrA(x,y,n)=(l + A) (l _ ^)-“ (^ + A)-(”+i-™) 

xiJ) (1 - re)™-' {w + A)”-"* ^ ) . 

\Jc„ 2z7ru; J 


Note that, by Gauchy’s residue formula, the second integral is equal to the coefficient 
of in (1 —t(;)'"“^(w-|-A)"“”’', which is a Krawtehouk polynomial. The first 

integral is of a similar nature, except that the role of the zeros {1, — A} and of the 
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poles {0, 00 } in the integrand are exchanged. This suggests to perform the change 


of variables 2 = 


i+A- 


(93) 


(1 + A)-"A" 


which transforms the first integral in (92) into: 

—du 


c„ 


‘(1 + A- 




2iTTU ’ 


where is a small negatively oriented contour encircling u = 0 (note that the 
image of the negatively oriented contour Cu under the change of variable z = 
is a negatively oriented contour encircling z = 1, which is homotopic in C U { 00 } \ 
{1, —A} to the positively oriented contour C^). Using Cauchy’s residue formula, the 
second integral is equal to the coefficient of in (1 + (1 — . 

Redistributing powers of A we finally obtain: 


Proposition 21. The biased creation rate reads 

(94) CrA(x, j/,n) = f j cx{A,B,n-l)cx{B,A,n-1) 

where A = n — m — y = {n — 1 — x — y)/2, B = m — 1 = {n — 1 + x — y)/2 and 
c\{A, B,n) is the coefficient of z^ in {1 — zffi{1 + z)^~^. 


Note that = (p — 1) and A ^ so we recover [CEP961 Proposition 23], 

up to the exchange p 1 —>■ 1 — p and the antidiagonal reflection {x + y,x — y) 
{—X — y,x — y) that correspond to the fact that we consider west-going rather than 
north-going dominos. |CEP961 Proposition 2] also follows, by taking A = 1, i.e. 
P = 1/2. 

We now turn to the edge-probability generating function. We need to multiply 
( [8^ by u^v'^t'^ = {tu)'^{t/uffi'^^~'^{uvYft and sum over all n > 0, to € [l..n] and 
p € Z. Assuming that A < 1, and can be taken as circles with center 0 and 
radiuses between A and 1. For t small enough and u,v of modulus close to 1, it is 
possible to interchange the double sum over n, to and the double contour integral 
to yield 


(95) 


Ux{u,v,t) = ^ 

y& 


1 

(2i7r)2 




(j){z, w)(uvw/zY 


dz dw 
z w 


where 


tu(l — 'w) t(l+X/w) 

(96) Mz w) = 

l-z u(l-|-A/z) ' ' 

The sum over all y allows to get rid of one integral. More precisely, assuming 
juuj < 1 and jraj fixed, we take two different contours Cz depending on the sign 
of p, namely a circle (resp. C~) of radius slightly larger (resp. smaller) than 
\uvw\ if p > 0 (resp. p < 0). Splitting the sum over p in two accordingly, we may 
interchange each sum with the integral, resulting in 


nA(u,u,t) = 


(97) 



1 / .dw 

= z—(7) (b{uvw,w) —. 
2*7’" /c„ w 


In the latter integral, the integrand has two poles, but only one of them falls within 
Cyj for t small and juuj close to 1, and we end up with: 


Proposition 22. The biased edge-probability generating function reads 
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For A = 1, we recover the expression given in |Dull[ IBPllj (again up to the 
change of variables needed to pass from west-going to north-going dominos). 


Remark 23. The generating function associated with the biased creation rate takes 
a much more symmetric form, namely 


(99) 


Clx{x,y,n)u^v'^t'^ = “^—(1 “ t/u)Ilxiu,v,t) 

x,y,n 

_ (A + l)t __ 

((1 -I- A)(l -I- - t{u + M“l) - Xt{v + ’ 


and it remains the same for other types of dominos. The combinatorial explanation 
of this phenomenon (and of the meaning of the term “creation rate”) comes from 
the domino shuffling algorithm [EKLP92bj, which implies that Gi\(x,y,n) is the 
probability that, in a biased random tiling of the Aztec diamond of size n, the 
2x2 square centered at (cc, y) is covered by exactly two dominos (regardless of 
their orientation). 


6.2. The inverse Kasteleyn matrix. Chhita and Young gave in |CY14[ Sec¬ 
tion 5] a formula for the inverse Kasteleyn matrix of the Aztec di amo nd of size n, 
for the biased distribution (whose definition is recalled on page 33). Under this 
distribution the probability for a tiling is proportional to q to the number of flips 
from the fundamental configuration times xfX to the number of vertical dominos. 
We now explain how to relate their formula with the entries 

Let a = a^) be an even vertex and /3 = {13^, 13^) an odd vertex, which have 

survived the contraction of edges. The coordinates we will use are those of the 
contracted graph, so G [0..2n — 2] and /3^ G [—1..2n — 1]. 

In this particular case, for k G [—1..2n — 1], the function Fk{z) from Equation ^ 
becomes: 


( 100 ) 


Fkiz) = 


1 


(1 - VXq-^^z) n;=L(fe+i)/ 2 j (1 + VXq^^+^/z) 


Chhita and Young use coordinates (xi, X 2 ) (resp. (j/i, 2 / 2 )) to localize odd (resp. even 
vertices. In their terminology, they are white and black respectively. These coordi¬ 
nates correspond to axes that are along the diagonals of Figure on the right. 

The correspondence between the two systems of coordinates is 


( 101 ) 


=1 + 13^ + 2/3y 
X 2 = 1 + ( 3 ^ 


yi = 1 -I- a’' -I- 2 q!^ 
2/2 = 1 - 1 - 0 ;’' 


Performing the change of variable 


( 102 ) 


c 


qV2-2 

w 


qy2-2 

z 


in the integral defining in [CY141 Theorem 5.1] in the case when X 2 < 2 / 2 ) 
one recovers the same factors for the rational fraction in ^ and 9 as for ^)) 

up to possibly numerical multiplicative constants. One has just to check that 
the contours enclose the same sets of poles. Under the change of variables, the 
contour Fi^^ becomes a contour for Q enclosing —'/Xq^^, —VXq ^^~^'^,..., —V~Xq'^^~^ 
(and which may or may not enclose zero, since the original rational fraction in w is 
regular at infinity), and the contour Fq is mapped to a large contour for 6 containing 
separating infinity from a domain containing all the poles and the contour for 9. 
This contour can be deformed freely as long as it does not cross the one for C, or 
infinity, because \/Fk(9) is a polynomial, and thus has no poles. 
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In the case when > 1 / 2 , the extra term in |CY14| comes from the residue at 
w = z, which can be integrated to the double integral, at the cost of interchanging 
the nesting of the contours. 


6.3. The arctic circle theorem. Under the uniform measure, domino tilings of 
a large Aztec diamond exhibit a spatial phase transition, known as the arctic circle 
phenomenon [,TPS9811.TohObj . Outside the inscribed circle, with probability expo¬ 
nentially close to 1, all dominos are arranged in a brickwall fashion. This is called 
the frozen region. Inside the circle, the probability of each orientation is non degen¬ 
erate and does not go to 0 or 1. We now explain how to recover this phenomenon 
from our formalism. The parameters q and A are now set to 1. 

Let a = = {2m,y — be an even vertex. The probability pa that a 

dimer connects a with {3 = {(3^,(3^) = (2m — is given by (89) with A = 1. 

In the scaling limit n —> 00 , m/n —> r, y/n —>■ y, this probability becomes 


(103) 


where 


Pa 


C, JC, 


exp(n(S'(z, r, x) - S{w, r, x) + o(l))) 


dzdtc 


{z — w){2iTry 


(104) S{z,T,x) = -rlog(l - z)-{l- 2 :)log(l -k 1/z) - xlog(z). 

We now proceed as in |OB03| to obtain the asymptotics of this probability. For 
fixed {t,x), the function z 1 —>■ S{z,t,x) has two critical points. 

• If the two critical points are real, the integral goes to 0 or 1 exponentially 
fast with n by the saddle point method. The point (r, x) is in the frozen 
region. 

• If the two critical points are complex conjugate, one can move the contours 
so that they cross transversally at the complex critical points to apply 
again the saddle point method. By doing so, we pick the contribution of 
the residue at z = w, which gives the main contribution of the integral, 
giving a result strictly between 0 and 1. 

The transition between those two regimes correspond to the value of (r, y) for which 
the two critical points merge. This happens when the discriminant of the numerator 
of dS{z,T,x)/dz is equal to zero. This gives 

(105) (2t- 1)2-4(r-f x)(l - X - "t) = 0, 


which under the change of variables 


(106) 


f u = 2t + x 

1 = X 


corresponds to the circle 2(u — 1)^ -|- 2v'^ = 1 inscribed in the limiting square of the 
Aztec diamond, given by |m — l|-|-|z)| < 1. 


7. Conclusion 


We have introduced the rail yard graph dimer model, and computed its partition 
function and correlation functions. We point out that it corresponds essentially to 
the most general Schur process with nonnegative transition probabilities, see the 
discussion in |Borll[ Sections 1 and 2]: we handle an arbitrary finite number of 
“a” and “/3” specializations, and any other specialization can be obtained by taking 
suitable limits (in particular, to get the Poissonized Plancherel measure, one shall 
consider the “dilute” limit of RYGs, namely take a sign sequence of the form 
an arbitrary LR sequence, and a constant weight sequence zin, then let n —>■ 00 ). 

Many directions can be explored from here. By applying the random generation 
algorithms of [BBB+14 , we may generate large RYG dimer configurations, which 
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allows to observe limit shape phenomena as in the cases of (skew) plane partitions 
[OR031IOR07) and of the Aztec diamond, discussed above. RYG seem to allow 
for an even larger variety of singular points on limit shapes, and of corresponding 
limiting processes, which are currently under investigation. 

The appearance of the rational edge-probability generating function (98) in the 
context of the Aztec diamond (when summing over diamonds of all sizes) raises 
the question whether such rationality phenomenon may subsist for other types of 
RYGs. A natural idea is to consider RYGs with periodic LR, sign and weight 
sequences. Preliminary research indicates the rationality phenomenon occurs only 
in another case, namely skew plane partitions of “staircase” shape. In other cases, 
we obtain an algebraic, but not rational, generating function (algebraicity being 
expected from the very nature of our computations). 

In this paper we have obtained the correlations for pure RYG dimer configura¬ 
tions (by computing vacuum-to-vacuum expectation values of fermionic operators). 
Other types of boundary conditions can be considered, as in |BGGI4) where the 
corresponding partition functions were computed (and the extension to RYGs is 
straightforward). However, adapting the computation of correlation functions done 
in the present paper is not so easy, since it requires an adaptation of Wick’s for¬ 
mula. For arbitrary but fixed boundary conditions, we know from general facts, 
namely the generalized Wick theorem |AZ13| or the Eynard-Mehta theorem |BR05| , 
that correlations will still be of determinantal form, however it is not clear how to 
compute explicitly the propagator/determinantal kernel. Such computation could 
be done by Petrov [Pet 14] for some lozenge tilings, and we are looking for other 
tractable cases. Also of interest is the case of free boundary conditions (that is, we 
sum over all possible boundary states). When only one of the boundaries of the 
RYG is free (corresponding to symmetric RYG dimer configurations), the correla¬ 
tions are known to be PfafHan |BR05| , and we have found an adaptation of Wick’s 
formula which would allow for a computation similar to that done in the present 
paper, bypassing the use of the (Pfaffian analogue of) Eynard-Mehta’s theorem. 
Details should appear in a subsequent publication, see also |DFR12l IVull3[ IPanl4j 
for related results. When the two boundaries of the RYG are free, the nature of 
the correlations is unknown, even though the partition function can be computed 
following the lines of |BGGI4j . We believe they should be the Pfaffian analogues 
of correlations for RYGs with periodic boundary conditions, related to the periodic 
Schur process of [Bor 07]. 

Finally, a tantalizing question is whether it is possible to consider “interact¬ 
ing” deformations of our dimer models. Besides the directions mentioned in the 
conclusion of [BCCI4| , let us mention that fermionic techniques have been recently 
used, together with methods from constructive field theory, to prove rigorous results 
about interacting dimer models, see e.g. |GMT14j and references therein. Another 
intriguing fact is that identities arising from z-measures (which are instances of 
Schur measures) have found applications in the context of quantum integrable sys¬ 
tems |KKM+11| . 
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Appendix A. Commutation of bosonic and fermionic operators 


In this section, we give a self-contained combinatorial proof of the commutation 
relations between the bosonic and the fermionic operators stated in Proposition [TT| 


Proof. Let us first prove (371. It is equivalent to prove that for any k and for any 

A, 


(107) 

Define 


rfl+(a::)'0fc|A) = {ipk + xipk-i)^R+{x)\X). 


(108) Uk = #{j > fc, Xj = •}. 

Case 1 : X has a white marble in position k. 

ri?,+ (a::)'0/c|A) enumerates the admissible dimer covers of an elementary rail yard 
graph of type (i?, -I-) and with right boundary equal to V'felA). Each dimer cover is 
specified by the value p. of the left boundary. Here, 

(109) i/>fe|A) = (-l)-MA(")), 

so the vertex in position k on the right boundary has to be incident to a certain 
edge Cfc in the dimer cover. Thus T]i.^.{x)'lpk\X) is a sum of two types of terms, the 
first (resp. second) type corresponds to horizontal (resp. diagonal). 

By the same argument as for the localization of horizontal dimers on a “double” 
column, dimer covers of the first type are in bijection with (and have the same 
weight as) dimer covers with right boundary A and with left boundary Such 
covers are enumerated by r/j+(x)|A). To obtain the original left boundary, /i, from 
this new cover, we need to apply ifk to rfl'+(a:)|A). The sign appearing is again 
(—!)”'“, because the number of black marbles above position k is the same in 
as in A. So the first term in the sum is equal to 

( 110 ) i;krR+ix)\X). 

By using the bijection used to localize diagonal dimers, and observing that the 
weights differ by a factor x, we obtain that the second term of the sum is equal to 

(111) Xlpk-l^ R+{x)\X). 


Note that the signs cancel out correctly because the number of black marbles 
above position fc — 1 on the left boundary is equal to the number of black marbles 
above position k on the right boundary. 

Thus we conclude in case 1. 

Case 2 : X has a black marble in position k. 


Here, i/'fclA) = 0, so the left-hand side of (1071 vanishes. 


rfl+(a;)|A) enumerates the admissible dimer covers of an elementary rail yard 
graph of type (i?, -I-) and with right boundary equal to A. Each dimer cover is 
specified by the value p, of the left boundary. The vertex in position k on the 
right boundary has to be incident to a certain edge in the dimer cover. Thus 
r/j_|_(x)'0fc|A) is a sum of two types of terms, the first (resp. second) type corre¬ 
sponds to Cfc horizontal (resp. diagonal). 

If /I. is a term of the first type, p has a black marble in position k, thus ifk + 
xtpk-i\p) = xipk-i\p) and 


( 112 ) xipk-i\p) = 




if p has a white marble in position k — 1 
otherwise. 


The sign is due to the fact that the number of black marbles in p above the 
position fc — 1 is equal to the number of black marbles in A above the position k, 
to which we must add the black marble in A in position k. 
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If is a term of the second type, /i has a black marble in position k — 1, thus 
i’k + = V’fclli) and 


(113) V'/cIm) = 


(—1)”'=/4(^') if fj, has a white marble in position k 


0 


otherwise. 


So all the nonzero terms of {ipk + xipk-i) rR+(a;)|A) have black marbles in posi¬ 
tions fc — 1 and k, and each term appears twice, with the same weight (because if 
/i is a term of the second type in r/j_|_(x)|A), it already carries a factor x coming 
from the diagonal dimer Cfc) and with opposite sign. 

Thus the right-hand side of (107) also vanishes. 


This concludes the proof of (37). Formula (41) is proved similarly, replacing 
ipk-i by ipk+i and zhy 

We noted that rfl+(x) was conjugated to ri+(a;) via w. Observe now that 4’* 
is conjugated to V'i via w, up to a sign (—1)'*' verifying (—= —(—1)®' (s; is 


defined by an equation analogous to (35)). This enables us to deduce formulas 
and 


Using that rij+(a;), Tlj^{x) and ijj{z) are respectively adjoint of rj{_(a;), rL_(a:) 
and ^/>(l/z), we deduce the last four formulas. 

□ 


Appendix B. Wick’s formula 
In this section, we provide a proof of the identity 


(114) (0|r (4'(ai), 'k*(/3i),..., 'k(a«), vI/*(^,)) |0) = 

det (0|r(vI/(a,),^'*(/^4)) 10) 

1<2J<S 


used in (59). 

Let n be the set of partitions of {1,..., 2s} into unordered pairs. An element 
T e n can be canonically written r = {{ii, Ji}, •. ■, {*si Js}} with A < *2 < • • • < *s 
and it < jt for all t, and viewed as a permutation 


(115) 


fl 2 ••• 2s-1 2 sA 

Vi h ■■■ is jsj ’ 


which allows to define its sign €{t). Recall that the Pfajfian of a 2s x 2s antisym¬ 
metric matrix A = (Aij)i<ij< 2 s is defined as 


S 

(116) Pf(Al) = ^ e(T) A(2t-i),r(2i)- 

T^n 


Let F be the space of (countably infinite) linear combinations of ipk’s and i/’^’s. 
Proposition 24 (Wick’s formula). For Xi,..., X 2 s elements of F, we have 
(117) (0|Ai.--A2«|0) =Pf(A) 

where A is the 2s x 2s antisymmetric matrix such that Aij = (0|AiAj|0) for i < j. 

Proof. Let F+ (resp. F~) be the vector space spanned by the ifk with fc > 0 and 
the tpl with fc < 0 (resp. the ifk with fc < 0 and the ipl with fc > 0), so that 
F = F+ 0 F~ . For X G F, we denote by A+ and X~ its projections on these two 
subspaces: observe that (0|A+ = 0 and A“|0) = 0. 
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Let Xi,..., X 2 s be elements of F. Observe that, by proposition {X-,X+}is 
a scalar for any i < j, and 


(118) 


m,x,\%) = (0|x+x 
= mx+x. 


■x+x- + X-X- + x-x+\%) 

■ x+x- + x-x- + {X-,X+} - x+x- 10) 


= {X-,X+}. 


We proceed similarly to compute the left-hand side of (117): we write each of the 
Xi as Xi = X^ + X~, expand the product and get a sum of 2^® terms. In each 
of these terms, we move all the X~ to the right. To get a nonzero contribution, 
we must pair up each X~ with one X'^ such that i < j, and multiply all the 
anticommutators obtained in this fashion. So each nonzero contribution is equal, 
up to a sign, to 


(119) 




t=i 


where r € 11 is some partition into pairs of {1,..., 2s}. Furthermore, the sign of 
this contribution corresponds to the number of swaps needed to bring each 
exactly to the left of i.e. it is the signature of r. This yields exactly the 

right-hand side of (116) with Aij = {X~,X^} = {^X^Xj\^). 


□ 


To prove the wanted identity ( |114| ), we set F 2 t-i = 4'(at) and T 2 t = 'l'*(,0t), 
t = 1,..., s. By the definition of natural ordering, there is a permutation a such 
that 


( 120 ) 

By Proposition 


24 


r(Ti,...,T2.) = e(a)l(,(i)---h;(2.). 

and since = T{Y^{i),Y„{j)) for any i < j, we have 


( 121 ) 


(0|r (Yi,..., r2s) | 0 ) = e(a) Pfi<,,,<2.(0|r(Y,(q, Y,(,))|0) 

= Pfl<z.,<2s(0|TO,Y,)|0) 

(when we simultaneously permute the rows and columns of an antisymmetric ma¬ 
trix, the PfafEan is multiplied by the sign of the permutation). Finally, we observe 
that {%\T(Yi, Yj)\%) is nonzero if and only if i,j have different parities, which allows 
to rewrite the Pfafhan as the wanted determinant without sign. 
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